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ABSTRACT 

The  relations  which  describe  the  mode  of  occurrence 
of  a  fresh  water-salt  water  system  in  a  coastal  aquifer 
are  reviewed  and  discussed.  For  comparison  two  models 
were  used  to  study  this  two  density  system  in  the 
laboratory,  and  from  the  studies  a  universal  set  of  curves 
were  derived,  which  relate  the  rate  of  flow  of  fresh  water 
into  the  aquifer,  the  permeability  of  the  pervious  material 
and  the  depth  to  the  interface  in  such  a  system.  Finally, 
as  a  discussion  of  the  laboratory  studies,  the  results  were 
applied  to  a  case  history  which  has  been  of  some  concern  to 


the  author. 
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c  c  = 

cubic  centimeters 

c  f  s  = 

cubic  feet  per  second 

e  = 

void  ratio 

G  = 

specific  gravity  of  soil  solids 

g.p -d.  = 

gallons  per  day 

H 

Distance  from  shoreline  of  in  flow  of  fresh  water 

Hi  = 

Distance  from  shoreline,  of  pumping  source 

K  = 

P  er  meab i 1  it  y 

L  = 

Length  of  island  portion  of  model 

min  .  = 

minute 

M  .  S  .  L  .  = 

Mean  Sea  Level 

ppm  = 

parts  per  million 

Q  = 

Quantity  of  flow 

Q?  = 

Quantity  of  flow 

q  = 

Quantity  of  flow  per  foot  width  of  aquifer 

x  = 

Coordinate  distance  from  shoreline 

x  = 

Coordinate  dist  an c  e 

Y  = 

Depth  to  interface 

y  = 

Coordinate  distance 

w 

unit  weight  of  water: 

II 

Specific  gravity  of  fresh  water 

p  - 

Specific  gravity  of  salt  water 
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INTRODUCT ION 

From  Biblical  times,  man  has  made  use  of  the 
tremendous  reserves  of  water  stored  in  the  ground.  This  has 
been  especially  true  in  countries  whose  geology  and 
topography  dictate  the  scarcity  of  surface  sources  of  water, 
such  as  lakes  and  rivers,  and  which  have  been  forced  to 
become  more,  dependent  on  ground  water  sources.  With  increased 
demands  for  water,  have  come  the  problems  of  more  accurate 
methods  of  estimating  ground  water  reserves.  For,  unlike 
surface  waters  where  the  storage  reservoirs  are  there  to  be 
seen,  and  simple  conventional  computational  methods  can  be 
used  to  estimate  the  volume  in  storage,  ground  water,  being 
stored  underground,  presents  much  more  of  a  problem  in  having 
its  storage  volume  estimated. 

More  recent  developments,  such  as  the  realization 
that  ground  water  sources  are  potentially  freer  from  con¬ 
tamination  of  radioactive  fallout,  have  led  to  a  greater 
interest  in  the  development  of  ground  water  storage  basins  as 
sources  of  water  for  public  supply. 

This  study  limits  itself  to  a  study  of  the  occurrence 
of  ground  water  in  an  aquifer  in  which  particular  geological 
conditions  manifest  themselves. 
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1 .  1  Ground  Water  in  _a  Coastal  Aquifer 

The  mode  of  occurrence  of  ground  water  in  certain 
coastal  areas  is  unique  and  of  special  interest.  The 
geological  cross-section  considered  in  this  study  shows  a 
relatively  permeable  stratum  such  as  sand,  limestone  or 
volcanic  tuff,  varying  in  thickness  from  a  few  feet  to 
hundreds  of  feet,  possibly  underlain  by  an  impervious 
stratum.  The  rainfall,  assumed  to  be  more  intense  away  from 
the  shoreline,  percolates  into  the  ground  under  the  action 
of  gravity  until  it  reaches  the  saturated  zone,  and  then 
flows  toward  the  sea.  Due  to  the  tortuous  nature  of  the 
passageways  between  the  pores  of  the  permeable  material  and 
the  small  hydraulic  gradient,  the  velocities  of  flow  are  very 
small.  The  static  water  level  in  the  saturated  zone  near  the 
shoreline  is  somewhat  above  mean  sea  level  and  is  usually 
higher  at  the  point  of  concentration  of  inflow  into  the  aqui¬ 
fer  than  at  the  coastline,  with  a  resulting  hydraulic 
gradient  to  the  sea. 

The  porous  nature  of  the  top  stratum  allows  also  the 
passage  of  sea  water,  which  being  heavier  than  the  fresh 
water  forms  a  hy  dr  o  dyn  am  ic  a  1 1.  y  balanced  system  with  the  fresh 
water.  In  this  system  the  fresh  water  floats  above  the  sea 
water,  the  two  waters  filling  the  pores  of  the  porous  media, 
with  the  sea  water  form  in g  the  bottom  and  outer  boundary  of 
the  underground  fresh  water  storage.  The  depth  to  the  inter¬ 
face  of  the  fresh  and  salt  water,  and  hence  the  volume  of 
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fresh  water  in  storage,  depends  primarily  on  the  height  of 
the  free  surface  of  the  fresh  water  above  mean  sea  level  and 
also  on  the  relative  densities  of  the  two  systems.  Various 
other  factors  determine  the  volume  of  fresh  water  in  storage 
in  such  a  system,  and  these  are  considered  in  some  detail  in 
the  present  work. 

The  coastal  areas  of  many  countries,  such  as  some 
of  the  seaboard  areas  of  the  coast  of  Holland  (Ref  1),  the 
coastal  areas  of  California  (Ref  2),  and  several  islands 
such  as  Hawaii  (Ref  3) ,  the  Bahamas  (Ref  4)  and  Barbados 
(Ref  5)  depend  on  this  coastal  ground  water  for  public  water 
supp ly . 

1 . 2  Purpose  and  Scope  o f  the  Present  Study 

This  unique  mode  of  occurrence  of  ground  water  in  a 
coastal  aquifer  has  been  of  considerable  interest  to  the 
author  while  carrying  out  field  investigations  in  Barbados. 
There  appeared  to  be  certain  discrepancies  in  the  values 
obtained  as  to  the  depth  of  the  interface  of  the  fresh  and 
salt  water  found  from  equations  describing  the  balance  of  the 
two  density  system  (Ghyb en-Her zb er g  relation).  Moreover  an 
extensive  literature  survey  on  studies  previously  made  on 
the  topic,  showed  that  no  quantitative  studies  were  under¬ 
taken  wh ich  related  the  inflow  rates  of  fresh  water  into 
the  aquifer,  the  position  of  the  interface,  the  permeability 
of  the  material  in  the  aquifer  and  any  other  parameters 
which  manifest  themselves  in  such  a  two  density  system. 
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This  study  attempts  by  means  of  models  of  a  fresh 
water-salt  water  system  to  provide  such  a  quantitative 
relation  and  in  general  to  study  the  limitations  of  the  use 
of  the  Gh yb en -H er zb er g  relation  in  both  a  pumped  and  a  non¬ 
pump  ed  system. 

Finally  as  a  discussion  of  the  results  obtained  in 
the  laboratory,  some  of  the  results  of  the  laboratory 
research  are  applied  to  an  actual  prototype  of  the  two 
density  system,  which  had  previously  been  studied  by  the 
author . 

1 . 3  Methods  o  f  Ap  proach 

The  work  has  been  divided  into  sections,  the  first 
of  which  briefly  discusses  the  history  of  the  Gh yb en -He r zb er g 
relation,  and  more  recent  additions  to  this  relation.  The 
first  section  also  deals  with  a  brief  discussion  of  studies 
on  salt  water  intrusion  into  such  a  system,  and  a  resume  of 
the  literature  on  some  of  the  studies  of  preventing  salt 
water  intrusion  is  given.  Two  brief  case  histories  are 
given,  one  dealing  with  sea  water  intrusion  along  the 
California  coast  and  the  other  dealing  with  the  occurrence 
of  ground  water  in  Barbados. 

The  second  section  of  the  work  deals  with  the  research 
conducted  in  the  laboratory.  Two  models  were  built,  one 
serving  as  a  prototype  of  the  other,  and  various  studies 
conducted.  This  section  has  been  further  subdivided  into  a 
quantitative  part,  and  a  qualitative  part  in  which  studies 
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were  conducted  to  observe  the  variation  of  the  position  of 
the  interface  with  pumping. 

The  final  section  deals  with  a  discussion  of  the 
results  of  the  laboratory  research. 


CHAPTER  II 


DEVELOPMENT  OF  RELATIONS  OF  BALANCE  OF 
A  TWO  DENSITY  SYSTEM 

During  the  last  century,  many  observations  have 
been  made  on  the  mode  of  occurrence  of  a  balanced  two 
density  system.  One  of  the  first  relations  which  describe 
the  mechanism  of  the  fresh  water-salt  water  balance  is  the 
Ghyb en-H er zb er g  relation.  The  relation  is  developed  and 
discussed  in  this  chapter. 

2 .  1  H is  t  or  y  o  f  Ghyben-Herzherg  Relat  ion 

The  first  enunciation  of  the  principle,  that  fresh 
ground  water  in  a  coastal  aquifer  floats  above  salt  water 
because  of  its  lower  density  was  contained  in  an  article  by 
Badon  Ghyben,  published  in  1888-1889  (Ref  1) .  The  discovery 
was  made  as  a  result  of  a  well  drilling  program  near 
A ms t  er  dam. 

In  a  publicat  ion  at  Munich  in  1901  (Ref  6)  Herzberg, 
from  a  study  of  water  level  fluctuations  and  salinity  in 
wells  for  the  water  supply  of  parts  of  the  North  Sea  coast, 
deduced  the  theory  that  fresh  water  floats  above  sea  water, 
because  of  the  difference  in  specific  gravity.  Herzberg 
deduced  a  ratio  of  1*37,  between  the  thickness  of  the  lens 
of  fresh  water  above  mean  sea  level,  and  the  depth  to  the 
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fresh  water-salt  water  interface  below  sea  level.  Herzberg 
observed  that  the  relationship  did  not  hold  exactly  in  all 
cases,  as  it  was  strongly  influenced  by  the  grain  size 
distribution  of  the  dune  sands,  but  that  it  was  approximately 
correct . 

Pennink  in  an  article  published  in  1905,  (Ref  7)  to 
a  large  extent  substantiated  the  findings  of  Ghyben  and 
Herzberg.  Pennink,  too,  worked  in  the  coastal  dune  area  near 
Amsterdam,  and  made  a  comprehensive  hydrologic  study  of  the 
area.  In  his  article,  he  showed  geological  cross-sections 
and  hydrological  cross-sections  showing  fresh-salt  water 
boundaries  and  hydraulic  grade  lines  in  the  various  aquifers. 
Pennink  in  the  study  added  to  the  usefulness  of  the  Ghyben- 
Herzberg  relation,  by  showing  how  salt  water  can  eventually 
be  drawn  into  a  well  even  though  the  bottom  of  the  screen  is 
above  the  original  fresh  water-salt  water  interface. 

In  the  same  year,  1905,  d'Antrimont  (Ref  8)  described 
an  experimental  study  using  sand  to  simulate  the  Belgian 
coast,  and  fresh  water  and  salt  water  of  proportionate 
density.  He  found  that  the  depth  to  the  zone  of  contact  of 
the  fresh  and  salt  water  was  proportional  to  the  height  of 
the  fresh  water  surface,  agreeing  with  his  own  field 
findings  and  those  of  other  investigators.  d'Antrimont  also 
found  that  the  zone  of  diffusion  was  very  small  and  sharp, 
and  that  a  p r o t r ub er anc e  formed  on  the  contact  zone  when  the 
upper  fresh  water  was  artificially  drained. 
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2  .  2  Prob  lems  Assoc  iat  ed  with  t.  h  e  Two  Dens  it  y  System 


Besides  the  usual  problems  associated  with  any  ground 
water  system,  there  are  some  special  problems  inherent  in 
a  two  density  ground  water  storage  system,  such  as  in  a 
coastal  aquifer.  The  main  problems  are: 

(a)  The  difficulty  of  estimating  reserves  in  storage. 

This  difficulty  presents  itself,  since  for  the 
system  of  fresh  water  floating  above  salt  water  to 
exist,  there  must  always  be  a  hydraulic  gradient  of 
the  top  fresh  water  surface  in  the  saturated  zone  to 
the  sea,  and  hence  fresh  water  is  being  lost  to  the 
sea  at  all  times.  Preliminary  investigations  have 
also  shown  that  near  the  shoreline,  there  is  a 
relatively  large  seepage  area,  extending  some  distance 
out  from  the  shoreline,  where  fresh  water  is 
continually  being  lost  to  the  sea  (Refer  Figure  40, 

App  end ix  B )  . 

(b)  The  danger  of  salt  water  intrusion.  The  porespaces 
of  the  top  permeable  stratum  act  as  conduits  for  the 
passage  of  fresh  water  from  precipitation  on  the 
land  and  also  as  conduits  for  sea  water  from  the 
adjacent  ocean.  In  cases  where  there  is  a  rapid 
withdrawal  of  fresh  water  by  pumping  from  the 
aquifer,  there  is  the  danger  of  salt  water  being  drawn 
into  the  pumping  well. 


2 . 3  Per  ivat  ion  o  f  Ghyben-Herzberg  R  e 1  at  ion 


FIGURE  1 

FRESH  WATER  IN  A  COASTAL  AQUIFER 


Figure  1  shows  an  idealized  section  through  a  small 
island.  Curve  BCE  represents  the  free  water  surface  of  the 
ground  water  in  the  saturated  zone  of  the  permeable  stratum, 
the  point  of  concentration  of  rainfall  being  assumed  to  be 
at  the  centre  of  the  island  G. 

Line  BDE  represents  mean  sea  level,  and  curve  BF A 
represents  the  interface  of  the  fresh  and  salt  water. 

h  =  CD  =  height  to  point  of  atmospheric  pressure 

in  fresh  water  existing  above  mean  sea  level 
at  point  being  considered. 

H  =  DA  =  depth  below  mean  sea  level  to  the  interface 


of  fresh  and  salt  water 


'•  - 


Ill 


P  =  Specific  gravity  of  fresh  water 

P  =  Specific  gravity  of  sea  water 

H  p  =  (H  +  h)  p 

p 

H.  =  •— l  °-  .  Equation  1 

h  ^ 

Hence  the  ratio  of  heights  of  fresh  water  existing 
above  sea  level,  and  the  depth  to  the  interface  below  sea 
level,  is  the  same  as  the  ratio  of  the  specific  gravity  of 
fresh  water  and  the  difference  in  specific  gravities  of  the 
salt  and  fresh  water. 

2 . 4  Assumpt ions  Made  in  Per  ivat  ion  of  Ghyben-Her zber g  R  e  1  at  ion 

Several  assumptions  have  been  made  in  the  above 
derivation  and  are  discussed. 

(a)  It  is  assumed  that  the  fresh-salt  water  system  is  in  a 
state  of  statical  equilibrium.  This  implies  that  no 
flow  is  taking  place,  but  it  is  obvious  that  this  is 
incorrect  since  from  density  considerations  only, 
without  flow  a  horizontal  interface  would  develop. 

Figure  40  in  Appendix  B,  indicated  the  boundary 
conditions  of  the  interface  at  the  shoreline  and  gives 
some  idea  of  the  way  fresh  water  is  lost  to  the  sea 
continuously.  Considering,  however,  the  very  slow 
movement  of  the  fresh  water  towards  the  sea  in  the 
order  of  one  foot  per  day  the  assumption  of  a  statical 
system  appears  to  be  reasonable. 

(b)  The  effect  of  the  diffusion  of  salt  water  into  the 
fresh  water  has  not  been  considered  in  the  derivation 


of  the  formula.  It  was  assumed  that  the  interface 


was  sharp.  Nomitsu,  Toyohara  and  Kamimoto  (Ref  9) 
have  shown  that  the  velocity  of  diffusion  between  the 
salt  water  and  fresh  water  is  very  small  compared  to 
the  horizontal  velocity  of  flow  of  fresh  water  above 
the  salt  water  and  hence  produces  only  very  small 
variations  in  the  interface  derived  by  the  use  of 
the  Ghyb en-Her zb er g  relation. 

2 . 5  Impr  ovement  on  Basic  Ghyben-Her zber g  R  e 1 at  ion 

In  a  paper  published  in  1927,  Nomitsu,  Toyohara  and 
Kamimoto  (Ref  9)  by  a  consideration  of  Darcy's  Law  of  flow 
through  porous  media  derived  equations  for  the  curves  BGC 
and  BFE  in  Figure  2. 

Equation  to  BGC  is 


y 


l2 


S* 


Equation  2 . 


and  equation  to  BFE  is 
y2  ~  H2 

R  X  .  Equation  3  . 

where  y,  y  ^ ,  x,  h,  H  and  R  are  defined  in  Figure  2. 

In  the  same  study  the  authors  derived  an  equation 

which  allowed  the  permeability  of  the  medium  through  which 

flow  is  taking  place  to  be  determined. 

Th  e  equat ion  is 

K/Q-csf  2R  P6 

b  H  2  ’  p~  ....Equation  4. 

Equations  2  and  3  indicate  the  parabolic  nature  for  the  curves 
BGC  and  BFE. 


To  apply  these  equations,  it  is  required  only  that  the 
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height  of  the  free  surface  of  the  fresh  water  in  the 
saturated  zone  above  mean  sea  level,  and  the  corresponding 
depth  to  the  interface  at  one  point,  as  well  as  the  horizontal 
distance  from  the  shoreline  at  the  point  being  cons ider ed,  b e 
known . 


FIGURE  2 

TOP  FRESH  WATER  SURFACE  AND  INTERFACE 
IN  A  COASTAL  AQUIFER 


Curve  BGC  shows  free  surface  of  fresh  water  in  the 
saturated. 

Curve  BFE  is  the  interface  of  the  fresh  water  and 


salt  water. 


and  y  are  respectively  the  ordinates  of  curves 
BGC  and  BFE  at  a  distance  x  from  the  shoreline. 

h  is  the  measured  height  of  fresh  water  above  mean 
sea  level  at  distance  R  from  shoreline. 

H  is  the  measured  depth  to  the  interface  from  mean 


sea  level  at  distance  R  from  the  shoreline. 


CHAPTER  III 


SEA  WATER  INTRUSION 

One  of  the  problems  encountered  when  the  ground  water 
stored  in  a  coastal  aquifer  is  used  as  a  source  of  public 
water  supply,  is  the  danger  of  polluting  the  reservoir  with 
intruded  sea  water,  if  the  safe  extraction  rate  of  the  fresh 
water  is  exceeded. 

Sea  water  intrusion  occurs  in  one  of  two  ways:  either 

as  a  wedge-shaped  body  of  salt  water  extending  itself  land¬ 
wards  beneath  the  fresh  water,  or  as  a  cone  of  salt  water 
extending  upwards  to  the  point  of  extraction  of  fresh  water. 
Here  too,  the  relative  densities  of  the  fresh  and  salt  water 
dictate  the  height  of  the  cone  of  salt  water  and  the  length 
of  the  intruded  wedge. 


3 . 1  Height  o  f  Cone  o  f  Sea  W  at er  in  an  As  sumed  Statical 
System 

Figure  3  is  an  illustration  of  dr  aw down  in  the  free 
water  surface  due  to  a  well  being  pumped  at  A,  with  a  sub¬ 
sequent  rise  in  sea  water  in  the  form  of  a  cone  with  apex  at 
C  as  shown.  Density  considerations  show  that  the  height  of 
the  cone  CD  is  relative  to  the  depth  of  drawdown  as  follows: 


h \ 

p 


CD 


P. 

P-P. 


AB 


E  qu  at  ion 


being  the  density  of  fresh  water 
the  density  of  sea  water. 
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FIGURE  3 


CONE  OF  SALT  WATER  BEING 
BRAWN  INTO  A  PUMPED  WELL 
IN  A  COASTAL  AQUIFER 


3 . 2  Length  o f  Intruded  Sea  Water  W edge 

Todd,  (Ref  10)  from  a  consideration  of  Darcy's  Law, 
derived  an  equation  for  the  length  of  an  intruded  sea  water 
wedge,  which  depends  mainly  on  the  flow  rate  of  fresh  water 
into  the  aquifer  and  the  permeability  of  the  material.  The 
equation  is 


q  =  %  /A-g,  Kb£ 

v  D  L  ....  Equation  6. 

■o 

q  is  the  seaward  flow  of  fresh  water  per  foot  of  ocean  front. 


< 
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yO  is  the  density  of  fresh  water. 

P  is  the  density  of  sea  water. 

L  is  the  lenth  of  the  wedge  of  sea  water 

b  is  the  depth  to  the  interface  at  the  furthest  point 

landward  of  the  wedge. 

K  is  the  coefficient  of  permeability  of  the  pervious 
s  t  r  at  urn. 

3 . 3  Prevent  ion  and  C  on  tr  o 1  o  f  Sea  Water  Intrus ion 

With  the  ever  present  danger  of  contamination  of  the 
reservoir  by  intruding  sea  water,  and  yet  with  a  desire  to 
obtain  maximum  extraction  of  fresh  water  from  this  source, 
especially  in  areas  which  depend  either  wholly  or  to  a 
considerable  extent  on  ground  water  in  a  coastal  aquifer  for 
their  water  supply,,  much  attention  has  been  given  to  methods 
of  counteracting  sea  water  intrusion. 

A  considerable  part  of  these  investigations  has  been 
carried  out  along  the  California  coast  and  inland  bays 
(Ref  2).  From  these  studies  most  investigators  have  agreed 
on  five  possible  methods  of  control  of  seawater  intrusion. 
Banks  and  Richter  (Ref  11)  list  these  as  follows: 

1)  Raising  of  ground  water  levels  above  sea  level  by 
reduction  or  rearrangement  of  pattern  of  pumping  draft 
or  both--Se award  hydraulic  gradient  must  be  maintained 
at  all  times. 

2)  Direct  recharge  with  fresh  water  of  overdrawn  aquifer 
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to  maintain  ground  water  levels  at  or  above  mean  sea 
level. 

3)  Maintenance  of  a  fresh  water  ridge,  above  sea  level  along 
the  coast.  This  is  similar  to  direct  recharge,  but.  is 
in  a  localised  area. 

4)  Construction  of  artifical  subsurface  barriers.  This 
can  be  accomplished  by  the  use  of  a  puddled  clay 
cutoff  wall,  sheet  piling,  injection  of  an  emulsified 
asphalt,  cement  grout,  bentonite,  silica  gel  or  calcium 
acrylate. 

5)  Development  of  a  pumping  trough  adjacent  to  the  coast. 

Some  of  these  methods  have  met  with  varying  degrees  of 
success  when  applied  in  the  field;  the  geological  conditions 
of  the  area,  as  well  as  the  avail ability  of  water  for  re¬ 
charge,  determining  the  measure  of  success. 

3 . 4  Sea  Water  Intrus ion  Along  California  Coast 

Sea  water  intrusion  along  the  California  coast  and 
inland  bays  has  received  much  attention  by  investigators  and 
considerable  literature  has  been  published.  Simpson, 

(Ref  12)  in  a  paper  published  in  1946,  described  the  extent 
and  movement  of  sea  water  intrusion  into  the  Salinas  Basin. 
The  ground  water  there  is  in  a  confined  aquifer  with 
pressure  of  180  feet.  As  of  October,  1945,  there  were 
approximately  6000  acres  contaminated  by  sea  water  extending 
to  nearly  two  miles  inland.  The  rate  of  encroachment  of  the 
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sea  water  was  600  feet  from  August  1944  to  August  1945. 

Simpson  in  his  investigations  found  that,  the  maximum  distance 
of  intrusion  was  to  the  most  inland  position  of  the  trough 
in  the  pressure  surface  under  conditions  of  heaviest  draft, 
and  was  estimated  to  include  9200  acres. 

Banks  and  B  o  ok man  (Ref  13)  in  a  report  to  the  State 
Water  Resources  Board  of  California  summarized  the  status 
of  sea  water  intrusion  in  California  and  possible  methods 
of  control  of  sea  water  intrusion.  These  methods  included 
recharge  through  injection  wells  by  the  Los  Angeles  County 
Flood  Control  District,  surface  spreading  and  subsurface 
barriers  by  grouting  and  by  puddled  clay  cutoff  walls.  The 
report  recommended  that  only  the  methods  of  pressure  ridges 
above  sea  level  and  construction  of  subsurface  barriers 
were  worthy  of  consideration  at  that  time. 

Laverty,  Jordan  and  van  der  Goot  (Ref  14)  carried  out 
tests  for  the  creation  of  fresh  water  barriers  to  prevent 
sea  water  intrusion  along  the  West  Coastal  Basin  and  summarized 
their  findings  in  a  report  in  1951.  The  tests  led  to  the 
following  conclusions: 

(a)  The  creating  of  a  fresh  water  ridge  along  the  coast  was 
the  roost  promising  solution  to  the  intrusion  problem. 

(b)  Bacterial  slimes  would  form  and  clog  aquifers  being 
recharged  unless  flow  was  sterilized.  Sterilization 
could  be  accomplished  by  dosing  recharge  water  with 
chlorine  at  an  initial  dosage  of  15  ppm. 


:  l 
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Chlorination  should  be  commenced  immediately  upon 
r  ech  ar  ge . 

(c)  Recharge  flow  displaced  saline  water  up  to  500  feet 
from  the  injection  well,  and  pumping  tests  indicated 
little  or  no  mixing  of  fresh  and  salt  waters. 

( d)  Recharge  well  carried  flow  of  between  1  and  2  cfs  for 
a  five  month  period. 

(e)  Pressure  elevations  created  by  recharge  were  approxi¬ 
mately  proportional  to  recharge  rate,  and  the  creation 
of  a  fresh  water  ridge  at  least  10  feet  in  height  was 
feasible  in  the  area. 


CHAPTER  IV 


GROUND  WATER  IN  BARBADOS 

The  island  of  Barbados  in  the  West  Indies  provides  a 
good  example  of  the  Ghyb en-Her zb er g  relation  in  ground  water 
contained  in  a  coastal  aquifer  and  its  inherent  danger  of 
sea  water  intrusion.  It  also  provides  a  good  basis  for  a 
discussion  of  the  results  of  the  investigations  carried  out 
in  the  laboratory  in  connection  with  this  thesis^  and  hence 
the  mode  of  the  occurrence  of  the  ground  water  in  Barbados  is 
discussed  here. 

4 . 1  Th e  Geo  logy  o f  the  Is  land 

Unlike  most  of  the  islands  of  the  Antilles  group  of 
the  West  Indies,  which  form  part  of  the  volcanic  arc 
associated  with  the  West  Coast  of  South  America,  Barbados 
and  Trinidad  belong  to  a  sedimentary  part  of  the  primary  arc 
extending  into  the  Brownson  Deep  (Ref.  15).  The  following 
summary  of  the  geology  of  Barbados  is  taken  from  S enn  (Ref. 20). 

During  the  Tertiary,  thick  s h a  1 1  ow-wa t er  elastics,  mud- 
volcano  deposits  and  deep-water  radiolarian  and  G  lob iger ina- 
rich  marls  with  volcanic  debris  accumulated  in  the  area  of 
Barbados.  This  sequence  was  interrupted  by  at  least  two 
orogenic  phases.  Orogeny  and  erosion  preceded  gentle  sub¬ 
sidence  during  the  Pleistocene  when  corals  became  established 
in  the  warm  shallow  waters  of  the  island  area.  Coral  growth 
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and  accumulation  of  sediment  continued  to  be  present  and 
varied  with  sea  level  fluctuations  due  to  glaciation  in  polar 
regions . 

Present  surface  geology  consists  of  porous  and  permeable 
coral  rock,  lying  on  impermeable  marls  and  shales  dipping 
predominantly  to  the  west  at  a  moderate  slope  (Figure  31, 
Appendix  B).  A  two  layered  hydrologic  system  is  thus  present, 
the  surface  topography  generally  following  the  topography  of 
the  underlying  impermeable  beds. 

4 . 2  Mode  o  f  Oc  cur r  enc  e  o  f  Gr  ound  Water 

Rainfall  which  occurs  predominantly  in  the  central 
higher  regions  of  the  island  averages  sixty  inches  per  year 
(Figure  32  in  Appendix  B).  It  is  estimated  that  one  fifth  of 
this  volume  of  rain  percolates  through  the  pervious  coral 
(Ref  16) ,  which  has  an  infiltration  rate  of  half  foot  per  day. 
As  the  water  moves  downward  under  the  action  of  gravity,  it 
eventually  reaches  the  impervious  clay,  where  its  downward 
movement  is  arrested,  and  the  water  than  travels  along  the 
dip  of  the  impervious  stratum  towards  the  sea.  Near  the  sea¬ 
shore,  the  heavier  sea  water  forms  a  fairly  effective  dam  to 
the  movement  of  the  fresh  water,  which  eventually  fills  all 
the  porespace  of  the  coral  to  a  height  somewhat  above  sea  level. 

The  height  above  sea  level  to  which  the  fresh  water 
rises  is  controlled  partly  by  the  permeability  of  the 
surrounding  coral,  and  partly  by  the  rate  of  inflow  of  fresh 
water  into  the  reservoir.  Finally  after  a  period  of  time 


, 
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the  system  reaches  the  hydrostatic  balance  of  the  Ghyben- 
Herzberg  relation.  This  coastal  reservoir  of  fresh  water 
has  been  given  the  name  of  sheetwater.  The  sheetwater  exists 
in  a  strip  along  the  entire  South  and  West  coasts  of  the 
island,  and  extends  landwards  to  the  point  where  the  clay 
stratum  rises  above  mean  sea  level.  The  sheetwater  basin 
varies  in  width  from  a  few  hundred  feet  to  about  four  miles 
in  places.  At  the  points  of  entry  of  the  flowing  ground 
water  into  the  coastal  system,  the  fresh  water  rises  to  a 
height  of  one  to  four  feet  above  mean  sea  level.  The  fresh 
water  as  it  moved  down  the  dip  of  the  clay  stratum  before 
reaching  the  sheetwater  system  described  above,  has  formed 
for  itself  over  a  period  of  years,  rather  discrete  channels 
by  solution  of  the  coral  along  the  coral-clay  interface. 

The  water  which  flows  in  these  discrete  channels  is  called 
"Stream  water"  in  Barbados.  Figure  33  in  Appendix  B  shows 
the  principal  directions  of  flow  of  underground  water  and 
the  approximate  demarcation  of  sheetwater  areas.  Figure  35 
in  Appendix  B  shows  a  section  through  the  island  showing  the 
-existence  of  the  ground  water  in  the  coastal  aquifer  as  it 
was  interpreted  by  S enn  (Ref  16).  Figure  39  (Appendix  B) 
was  prepared  from  field  measurements  taken  in  Barbados  by  the 
author  in  1961,  and  shows  the  correct  form  of  the  mode  of 
occurrence  of  the  two  density  ground  water  system  in  Barbados. 

4.3  Extraction  Sources 


Most  of  the  island's  population  depends  on  the 


government  water  supply  system  for  its  water.  Up  until  about 
1948,  only  the  stream  water  was  tapped.  Subsequently  as  water 
demand  increased,  and  the  existence  of  the  sheetwater 
determined,  three  pumping  sources  in  the  sheetwater  have  been 
developed.  Figure  34  in  Appendix  B  shows  the  location  of  the 
main  pumping  sources  for  the  island. 


Two  wells  in  the  sheetwater  area,  Belle  and  Haymans 


are  of  particular  interest  to  this  study.  The  Belle  well 
supplies  eight  million  gallons  a  day  out  of  a  total  of  twelve 
million  gallons  a  day  for  the  island's  water  demand.  Until 
April  1961,  Haymans  well  produced  two  million  gallons  a  day, 
but  this  caused  sea  water  intrusion  into  the  well  and  the 
yield  from  this  well  was  reduced  to  three-quarter  million 
gallons  a  day.  Haymans  well  possesses  a  graphic  history  of 
salt  water  intrusion  and  recent  observations  are  given.  These 
figures  are  taken  from  an  interim  report  by  Martin  Kaye 
(Ref  17)  to  the  Government  of  Barbados. 

Date  R  emar ks 


13/3/6  1 


Pumping  at  1.6  x  10^  g.p.d.; 

Water  level  in  well  4.8  inches  above  mean 
sea  level; 

132  ppm  chlorides  as  Sodium  Chloride  at 
pump  intake  which  is  two  and  one-half 
feet  below  surface  of  water; 

3550  ppm  chlorides  at  seven  feet  below 
water  surface  (bottom  of  well) ; 

At  this  stage  pumping  was  reduced  to 
850,000  g.p.d. 


14/3/61 


12  hours  after  reduction  of  pumping; 
Water  level  7.2  inches  above  mean  sea 
level ; 

94  ppm  chlorides  at  top  water  level; 


168  ppm  chlorides  at  two  and  one-half 
feet  below  surface  of  water; 

212  ppm  chlorides  at  five  feet  below 
surface  of  water; 

1050  ppm  chlorides  at  seven  feet  below 
surface  of  water  (bottom  of  well) . 
Pumping  at  1  x  1 0^  g.p.d.; 

Water  level  8.4  inches  above  mean  sea 
level ; 

88  ppm  chlorides  at  suction  of  pump; 

28/3/61 

Pumping  at  1  x  10 6  g.p.d.; 

Water  level  7.2  inches  above  mean  sea 
level  ; 

96  ppm  chlorides  at  suction. 

11/4/61 

Pumping  at  1  x  10^  g.p.d.; 

Water  level  7.2  inches  above  mean  sea 
level  ; 

68  ppm  chlorides  at  top  water  level; 

80  ppm  chlorides  at  suction  of  pu  mp  ; 
216  ppm  chlorides  at  four  and  one-half 
feet  below  the  surface  of  water; 

440  ppm  chlorides  at  seven  feet  below 
surface  (bottom  of  well). 

The  figures 

demonstrate  a  case  of  salt  water  intrusion. 

which  was  corrected  by  reducing  the  pumping  draft. 


CHAPTER  V 


LABORATORY  STUDIES  OF  A  FRESH  WATER-SALT  WATER  SYSTEM 

In  the  laboratory,  studies  were  carried  out  in  an 
attempt  to  reproduce  the  Ghyben-Her zberg  relation  and  to 
determine  its  applicability  under  vary  in g  inflow  rates  of 
fresh  water  into  the  system,  and  at  the  same  time  to  study 
the  effects  of  varying  rates  of  extraction  of  fresh  water 
from  the  system.  For  these  studies  two  models  were  used, 
in  which  the  permeable  strata  consisted  of  different  types 
of  s  and . 

5 .  1  Construct  ion  of  Mo  dels 

The  smaller  model  consisted  of  a  watertight  rectangular 
plywood  box  two  feet  high,  three  feet  long  and  three  and  one 
half  inches  wide  with  a  clear  transparent  plexiglas  front, 
one-quarter  inch  thick.  A  two  inch  square  opening,  placed 
eighteen  inches  from  the  bottom  in  one  of  the  small  sides  of 
the  box,  by  serving  as  a  fixed  outflow  weir  allowed  the  levels 
of  the  sea  water  and  of  the  fresh  water  which  overflowed  into 
the  sea  chamber  to  be  kept  constant.  Fine  standard  Ottawa 
sand  was  used  as  the  permeable  stratum,  while  the  bottom 
of  the  box  represented  the  impervious  stratum.  Figure  38  in 
Appendix  B  is  an  illustration  of  the  small  model. 

The  larger  model  was  similar  in  construction  to  the 
small  model,  but  with  dimensions,  sixteen  feet  long,  four  feet 


high  and  six  inches  wide.  The  fixed  outflow  weir  was 
placed  thirty-eight  inches  from  the  bottom  of  the  model. 

Fine  gopher  sand  was  used  as  the  permeable  stratum.  Figure 
37  in  Appendix  B  is  an  illustration  of  the  large  model. 

In  both  models,  the  sand  occupied  s  e  v  en  -  e  igh  t  h  s  of 
the  vo lume  of  the  box,  while  the  r  ema ining  one-eighth  was 
used  as  a  salt  water  chamber.  Fine  wire  screens  were  used 
to  prevent  the  sand  from  flowing  into  the  sea  chamber,  and 
to  retain  the  s  a  me  shoreline  conditions  throughout  the 
testing  program.  Clean  water  containing  enough  sodium 
chloride  to  give  it  a  specific  gravity  of  1.040,  and 
coloured  with  Fluorescein  dye  was  used  to  simulate  the  sea 
water.  The  salt  water  flowed  into  the  sea  chamber  from  a 
large  over-head  tank,  and  the  salt  water  in  the  sea  chambers 
was  maintained  at  a  fixed  level  and  fixed  density  by  a 
control  of  flow  over  the  weir  at  the  sides  of  the  models. 

The  specific  gravity  of  the  salt  water  was  determined  from 
hydrometer  readings.  Clean  uncoloured  fresh  water  was  dis¬ 
charged  on  to  a  specified  point  in  the  sand  from  a  constant 
head  over-head  tank.  Various  discharges  of  fresh  water, 
simulating  various  inflow  rates  into  the  aquifer,  were 
controlled  by  a  valve.  The  flow  rates  were  measured  by 
timing  the  flow  into  a  graduated  cylinder. 

To  simulate  pumping  conditions  in  the  field,  one-half 
inch  diameter  screens  of  fine  wire  mesh  were  inserted  into 
the  sand,  and  water  extracted  by  siphoning  through  rubber 
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tubing  at  a  controlled  rate.  A  small  outflow  of  salt  water 
was  maintained  over  the  weir,  and  this  water  surface  of  the 
salt  water  in  the  sea  chamber  was  used  as  a  datum  line. 

Level  lines  were  put  in  with  a  carpenter's  level.  A  vertical 
line  through  the  point  of  intersection  of  the  near  edge  of 
the  sand  and  the  free  water  in  the  sea  chambers  served  as  a 
reference  line  for  horizontal  measurements.  A  one  inch  and 
a  four  inch  square  grid  were  drawn  on  the  plexiglas  fronts 
of  the  small  and  large  models  respectively,  with  the  two 
reference  lines  discussed  above  as  axes  of  the  coordinate 
system.  All  measurements  of  the  heights  of  the  free  fresh 
water  surface,  the  distance  from  the  shoreline,  and  the 
depths  of  the  interface  were  referred  to  this  grid. 

5 . 2  P  hys ic  a  1  Properties  o  f  Sands  Used  in  Models 

C.190  standard  Ottawa  sand  was  used  as  the  permeable 
medium  in  the  small  model,  while  gopher  silica  sand  was  used 
in  the  big  model.  In  Appendix  C  are  given  the  results  of 
specific  gravity,  sieve  analysis  and  permeability  tests  which 
were  carried  out  on  the  two  sands.  The  tests  followed  the 
procedure  as  shown  in  the  manual  by  Lamb  (Ref  18) .  The 
specific  gravity  of  the  Ottawa  sand  was  found  to  be  2.67 
while  that  of  the  gopher  sand  was  2.68. 

A  sieve  analysis  showed  that  both  sands  were  rather 
uniform,  the  Ottawa  sand  having  97.4  per  cent  by  weight 
between  the  Number  16  and  Number  30  sieve  size  range.  The 
gopher  sand  had  68.7  per  cent  by  weight  between  the  Number  50 


and  Number  100  sieve,  while  29.1  per  cent  was  between  the 
Number  100  and  Number  200  sieve  size.  In  general,  then,  the 
mean  diameter  of  the  grains  of  the  gopher  sand  was  one 
standard  sieve  size  smaller  than  the  mean  diameter  of  the 
grains  of  Ottawa  sand. 

The  results  of  the  permeameter  tests  are  discussed  in 
detail  in  Section  6.7. 

5 . 3  T  es  t  ing  Pro  c  edur  e 

The  testing  program  was  divided  into  two  sections. 

One  section  was  devoted  to  detailed  quantitative  studies  of 
the  relation  between  the  position  of  the  interface,  the  flow 
rate  of  fresh  water  into  the  model,  the  permeability  of  the 
model,  and  the  distance  from  the  shoreline  to  the  point  of 
concentration  of  fresh  water  into  the  model. 

The  second  section  of  the  testing  program  was  more 
general,  and  attempted  to  show  changes  in  the  interface  under 
simulated  conditions  of  pumping.  The  distances  from  the 
shoreline  of  the  extraction  sources,  the  number  of  extraction 
sources,  and  the  extraction  rates  were  varied. 

In  order  to  conduct  a  test,  the  salt  water  of  correct 
density  was  made  up  by  the  addition  of  sodium  chloride  to 
fresh  water  contained  in  a  large  plastic  drum.  The  density 
was  checked  by  means  of  a  hydrometer.  Fluorescein  dye  was 
added  to  the  salt  water,  and  the  mixture  pumped  into  the  over¬ 
head  tank,  from  where  it  flowed,  its  rate  controlled  by  a  valve, 
into  the  sea  chamber  of  the  model.  The  level,  of  water  in 


the  sea  chamber  was  kept  constant  by  making  sure  that  a  small 
overflow  of  salt  water  from  the  sea  chamber  over  the  weir  in 
the  side  of  the  model  was  maintained.  The  salt  water  was 
collected  in  a  waste  tank  for  re-use.  At  the  same  time,  the 
required  flow  of  fresh  water  was  allowed  to  flow  on  to  the 
sand  of  the  model  at  a  predetermined  distance  from  the  shore¬ 
line.  The  rate  of  flow  of  fresh  water  was  determined  and 
checked  periodically  by  timing  the  flow  into  a.  graduated 
cylinder.  Hydrometer  readings  were  continually  being  taken 
of  the  water  in  the  sea  chamber ,  and  if  the  density  of  the 
water  tended  to  decrease,  due  to  any  mixing  of  fresh  water 
flowing  from  the  island  part  of  the  model,  then  a  larger  flow 
of  sa.lt  water  was  allowed  into  the  sea  chamber  and  hence  the 
density  of  the  salt  water  in  the  sea  chamber  was  kept  constant. 

Temperature  readings  of  the  fresh  and  salt  water  were 
taken  regularly,  but  throughout  the  whole  testing  program 
there  was  less  than  a  degree  change  in  temperature  of  the 
waters.  The  fresh  water  maintained  a  temperature  of  seven¬ 
teen  degrees  Centigrade,  while  that  of  the  salt  water  was  a 
constant  twenty  degrees  Centigrade. 

Soon  after  the  start  of  the  test,  an  interface  was 
clearly  seen  between  the  fresh  and  salt  water,  which  by  this 
time  had  both  percolated  into  the  sand.  In  the  case  of  the 
big  model,  two  days  were  allowed  for  each  run,  and  for  the 
small  model  one- half  day  sufficed  for  equilibrium  conditions 
to  establish  themselves,  that  is,  the  interface  to  remain 


fixed  in  a  certain  position  and  shape  for  that  particular 
flow  of  fresh  water  into  the  model.  This  position  of  the 
interface  was  recorded  by  reference  to  the  grid  which  was 
drawn  on  the  plexiglas  front  of  either  model.  A  final 
check  was  made  on  the  flow  rate  of  fresh  water  into  the 
model  and  this  was  recorded,  as  well  as  the  distance  from 
the  shoreline  of  the  point  of  application  of  the  fresh  water. 

In  order  to  determine  the  height  of  the  water  table 
above  the  level  of  water  in  the  sea  chamber,  one-sixteenth 
inch  diameter  tubes  were  inserted  into  the  sand,  thus 
creating  small  holes  in  which  the  free  surface  of  the  fresh 
water  was  clearly  visible  and  heights  above  sea  level  easily 
me  as  ur  e  d  . 

Several  runs  were  conducted  for  varying  inflow  rates 
of  fresh  water.  For  the  first  section  of  the  experimental 
work  in  which  quantitative  analyses  were  to  be  carried  out 
on  the  data,  the  point  of  application  of  the  fresh  water 
was  varied  twice.  In  one  case  the  ratio  of  the  distance  from 
the  shoreline  to  the  point  of  application  of  fresh  water,  and 
the  length  at  sea  level  of  the  island  portion  of  the  model, 
was  0.93  for  both  models,  and  in  another  series  of  tests  the 
ratio  was  0.43  in  the  small  model  and  0.54  in  the  big  model. 
For  the  second  part  of  the  testing  program  in  which  only 
general  qualitative  data  were  required,  the  experimental 
procedure  was  very  similar,  but  with  time  and  the  rate  of 
extraction  of  fresh  water  from  a  point  source  in  the  sand  of 
the  model  added  as  parameters. 
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After  equilibrium  conditions  had  been  obtained  for 
a  certain  flow  rate  of  fresh  water  into  the  model  similar  to 
the  tests  in  the  first  section  of  the  experimental  work,  a 
one-half  inch  diameter  cylinder  eight  inches  long  and  made 
of  fine  wire  screen  was  inserted  into  the  sand  to  a  depth 
two  inches  below  sea  level.  Thin  rubber  tub  in g  was  inserted 
into  the  well  thus  formed,  and  fresh  water  was  extracted 
through  the  tube  by  siphoning.  The  original  flow  of  fresh 
water  into  the  model  was  maintained.  Control  of  extraction 
rate  of  fresh  water  was  obtained  by  the  manipulation  of  a 
valve  on  the  siphoning  line.  The  rate  at  which  the  fresh 
water  was  being  extracted  was  measured  by  timing  the  flow 
into  a  graduated  cylinder.  The  time  at  which  extraction 
started  was  recorded,  and  the  different  positions  and 
shapes  of  the  interface  recorded  at  measured  time  intervals. 
The  rates  of  extraction  of  fresh  water  were  varied,  as  well 
as  the  point  of  extraction.  In  some  cases,  it  was  observed 
that  the  interface  developed  a  cone  shape  near  the  point  of 
extraction,  which  extended  itself  up  to  the  well,  and  that 
salt  water  was  being  extracted.  The  time  which  had  elapsed 
from  the  start  of  pumping  to  the  time  when  this  cone 
developed  was  recorded.  A  record  was  also  kept  of  the 
position  and  shape  of  the  top  free  surface  of  the  fresh  water 
as  they  changed  with  time  due  to  effect  of  extraction  of 
fresh  water.  In  some  of  the  tests,  the  position  and  shape 
of  the  interface  were  traced  on  the  plexiglas  front  at  hourly 
or  half  hourly  intervals  and  photographed. 


CHAPTER  VI 


RESULTS  AND  TREATMENT  OF  DATA 

In  this  chapter  the  results  of  several  trials  which 
were  conducted  in  the  two  models  are  given.  The  data  are 
arranged  in  dimensionless  groups  for  consideration.  The 
results  of  different  methods  of  determining  the  permeabilities 
of  the  two  sands  are  discussed. 

6 .  1  S  ymb  o 1 s  and  Terms  Used  in  P  r  e  s  en  t  at  ion  of  Data 

The  following  parameters  were  used  in  the  model  studies 
and  are  def  ined : 

p  and  p  are  the  specific  gravities  of  the  fresh  and  salt 
salt  water  respectively.  These  were  constant  at.  1.0  and 
1.043  respectively  through  the  majority  of  the  tests. 

L  is  the  horizontal  length  in  inches  of  the  land  portion 
of  the  model  measured  at  the  sea  level  line.  The  sea 
level  line  or  mean  sea  level  is  the  level  of  the  salt 
water  in  the  sea  chamber,  when  a  small  flow  over  the  weir 
is  taking  place. 

H  is  the  horizontal  distance  in  inches  of  the  point  of 
application  of  fresh  water  on  to  the  model  measured 
from  the  shoreline. 

X  is  the  horizontal  distance  in  inches  to  the  point  being 
considered  on  the  land  surface  from  the  shoreline. 

Y  is  the  vertical  distance  from  the  sea  level  line  in 


inches  to  the  interface  of  the  fresh  and  salt  water  for 


a  point  X  inches  from  the  shoreline. 

Yl  is  the  vertical  distance  in  inches  from  the  sea  level 
line,  to  the  top  free  surface  of  the  fresh  water  for  a 
point  X  inches  from  the  shoreline. 

Q  is  the  total  flow  rate  of  fresh  water  into  the  model  and 
is  measured  in  cubic  centimeters  per  minute. 

q  is  the  flow  rate  of  fresh  water  into  the  model  per  unit 
width  of  the  model,  and  is  expressed  in  cubic  centimeters 
per  second  per  centimeter. 

K  is  the  permeability  of  the  sand  in  each  model  and  is 
measured  in  centimeters  per  second. 

6 . 2  T  r  e  at  men  t  o  f  Data 

Tables  V  to  XXII  in  Appendix  A  show  the  values  of 
X,  Y,  Y 1 ,  and  H  for  the  different  rates  of  inflow  of  fresh 
water  into  the  models  for  the  first  section  of  the  experi¬ 
mental  work,  while  Figures  4  to  7  in  Appendix  B  show  the 
shapes  and  positions  of  the  interface  obtained. 

In  order  to  make  the  data  more  universal  in  its 
application,  dimensionless  parameters  were  used  in  presentation 
of  the  data.  This  also  made  possible  a  quick  and  easy  check 
on  the  applicability  of  the  results  of  the  laboratory  work, 
by  considering  the  small  model  as  truly  a  model  and  the  larger 
model  as  a  prototype.  This  was  desirous  since  no  reliable 
and  usable  data,  such  as  permeability  and  depths  to  inter¬ 
face  of  a  field  prototype,  were  available. 
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6 . 3  D Imens ion  less  Analys is 


The  parameters  ,  p  ,  H ,  L ,  K ,  q ,  X  and  Y  discussed 
earlier,  were  placed  in  dimensionless  groups  f)/p,  X/H,  Y/H, 
q/KH  and  H/L.  PJ  P  was  kept  constant  throughout  the  trials, 
so  that  X/H,  Y/H,  q/KH  and  H/L  were  considered  as  the 
variable  parameters.  q/KH  can  thus  be  considered  as  a  function 
of  X/H,  Y/H  and  H/L.  This  particular  choice  of  dimensionless 
grouping  appeared  to  be  the  most  favourable  one,  since 
inspection  of  the  data  obtained  from  the  tests  had  indicated 
the  importance  of  the  parameter  H.  For  future  considerations, 
another  dimensionless  group  q  T/YL,  where  T  is  time,  was 
established. 

6 . 4  Further  Tr  eatment  of  the  Data 

Preliminary  plots  on  arithmetic  paper  of  the  coordi¬ 
nates  of  the  interface  for  various  runs  indicated  a  parabolic 
shape  of  the  interface  (Figures  4  to  7).  This  parabolic 
relation  was  further  indicated  by  a  plot  on  double  logarithmic 
paper  of  X/H  versus  Y/H  for  different  values  of  q/KH .  These 
plots  produced  straight  lines  with  slopes  approximately  one- 
half.  However  it  was  observed  that  as  the  inflow  rate 
increased  in  either  model,  the  straight  line  which  resulted 
from  the  logarithmic  plot  of  the  data  in  the  form  X/H  versus 
Y/H  deviated  more  and  more  from  the  slope  of  one-half,  and 
that  these  tests  yielded  data  in  which  there  was  an 
appreciable  depth  to  the  interface  at  the  shoreline.  On  the 
arithmetic  plot  of  the  data,  it  was  easily  seen  that  with  the 


position  of  the  shoreline  as  the  origin,  most  curves  inter¬ 
sected  the  X  axis  some  distance  C  from  the  origin.  A  series 
of  parabolas  of  the  form  Y2  =  aX  were  drawn  on  transparent 
paper  for  a  wide  range  of  values  of  a  and  the  X  values 
within  the  range  observed  in  the  models,  and  these  curves 
fitted  to  the  arithmetically  plotted  data,  with  the  X  axis 
coinciding.  The  horizontal  displacement  of  the  Y  axis  of 
the  best  fitting  theoretical  parabola  from  the  Y  axis  of  the 
arithmetically  plotted  curve  obtained  from  experimental 
data  was  determined.  Thus  an  accurate  value  of  the  distance  C 
was  determined.  This  distance  C  was  then  algebraically 
added  to  the  X  ordinates  of  the  experimental  data  of  the 
curve  concerned.  This  in  effect  amounted  to  a  transformation 
of  axes,  giving  each  curve  a  different  origin. 

6 . 5  Pet  er minat ion  o  f  P  er meab i 1  it  y  K  From  Exp  er iment  a  1 
Data 

Equation  4,  section  2.5,  provided  a  method  for 
determining  the  permeability  K  of  the  sand  in  each  model 
for  different  inflow  rates  of  fresh  water.  In  the  equation 

K/Q  ££  2 X  Po  .  Equation  4. 

bY2‘  p_P' 

H0  is  the  value  of  Y  for  a  station  at  which  X  is  equal  to  R. 
The  equation  can  thus  be  more  generally  expressed  as 


K/Q  ca  _2X  Po  .  Equation  7. 

bY2’  P  -  P 
1  ‘o 

Thus  K  ^  2Q  P  .  Equation  8. 

bWx  ‘  p_po 


For  each  test  the  values  of  Y  2  /  X  were  determined  for 
several  values  of  X  and  corresponding  Y.  Values  of  X  and  Y 
for  points  near  the  shoreline  and  landwards  of  the  point  of 
inflow  of  the  fresh  water  were  not  considered,  for  reasons 
discussed  in  section  7.5.  An  average  value  of  Y2/X  was 
found  and  K  determined  by  substitution  of  the  appropriate 
values  into  Equation  8.  Finally  the  value  of  q/KH  for  each 
test  was  determined. 

Tables  XXIII  to  XL  in  Appendix  A,  show  the  X  and  Y 
values  of  the  data  of  Tables  V  to  XXII  referred  to  the  new 
axes,  C  being  the  distance  the  Y  axis  was  transformed.  The 
Tables  also  show  the  values  of  X/H,  Y/H,  Y2/X,  K  and  q/KH. 

6 . 6  Determinat ion  o  f  P  ermeab i 1  it  y  from  Velocity 
Measur ement s 

Two  tests  were  carried  out  for  flow  rates  of  fresh 
water  into  the  models  of  14  cubic  centimeters  per  minute, 
and  16  cubic  centimeters  per  minute  in  the  big  and  small 
models  respectively.  The  velocities  of  flow  of  the  fresh 
water  through  the  pores  of  the  sand  and  the  streamline 
pattern  were  determined  by  injecting  dye  at  various  depths. 
Figures  8,  9  in  Appendix  B  show  the  stre ami ines  obtained  in 

the  big  model  and  small  models  respectively  for  the  inflow 
rates  mentioned.  The  dye  was  injected  after  equilibrium 
conditions  had  been  established.  Tables  I,  II  give  the 
values  of  velocities  in  the  section  where  horizontal  flow  was 
evident,  and  also  the  values  of  permeability  obtained  in  the 


two  mode  Is . 


TABLE  I 


VELOCITIES  OF  FLOW  AT  DIFFERENT  DEPTHS 
IN  BIG  MODEL 


Depth 

below  sea 
in  inches 

level  Velocity 

in  feet  per  second 

Sea  leve 1 

1.35 

x  1 0  “  5 

7 

1.41 

x  10"5 

9 

1 .  1 

x  10" 5 

16  .  8 

2 .05 

X 

t— 1 

o 

i 

Ui 

Average 
K  = 

velocity  = 
7.02  x  10 ” 

1.5  x  10  feet  per  second. 

2  centimeters  per  second. 

T  AB  LE 

II 

VELOCITIES  OF  FLOW  AT 
IN  SMALL 

DIFFERENT 

MODEL 

DEPTHS 

Depth  below  sea  level 
in  inches 

in 

Veloc  ity 
feet  per  second 

2 . 5 

1.42 

X 

10-4 

4 . 6 

1 . 42 

X 

10"4 

7 

1.47 

X 

-4 

10 

Average  velocity  =  1.44  x  10-4  feet  per  second. 

K  =  .875  centimeters  per  second. 
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The  slope  or  hydraulic  gradient  of  the  interface  in 
a  short  section  where  the  streamlines  were  horizontal  was 
determined  for  each  model.  Since  the  e qu ip o t ent  ia 1  lines  in 
this  section  would  be  vertical,  it  can  be  assumed  that  the 

P 

slope  of  the  interface  divided  by  _ 2 _  gives  the  hydraulic 

P-P. 

gradient  of  the  top  fresh  water  surface  with  in  the  s  ate 
section.  From  Darcy's  relation  (Ref  10)  that  the  permeability 
is  the  velocity  divided  by  the  hydraulic  gradient  of  the  top 
water  surface,  and  with  the  assumption  that  the  measurements 
obtained  by  the  use  of  the  tracer  dye  gave  very  close 
values  to  the  actual  velocities  of  flow  in  the  models,  the 
permeabilities  were  calculated  to  be  6.0  x  10  "  2  cent ime  t  er  s 
per  second  and  0.875  centimeters  per  second  in  the  big  and 
small  models  respectively. 

6 . 7  Determinat ion  of  P  ermeab i 1 ity  From  P  er meamet er  Tests 
Following  the  procedure  described  by  Lamb  (Ref  18), 
permeameter  tests  were  carried  out  on  samples  of  the  sand 
used  in  each  model.  For  the  gopher  sand,  used  in  the  big 
model,  a  permeability  of  2.68  x  10“2  centimeters  per  second 
was  obtained  at  a  void  ratio  of  0.48,  while  for  the  Ottawa 
sand  in  the  small  model,  a  permeability  of  1.68  x  10 “2 
centimeters  per  second  was  obtained  for  a  void  ratio  of 
0.52.  The  data  for  the  tests  are  given  in  Appendix  G. 


6  .  8  Results  o f  Tests  in  Second  Sect  ion  o f  t_h e  Laboratory 


Studies 

In  addition  to  the  symbols  and  terms  defined  in 
section  6.1,  the  following  symbols  were  used  and  are  defined. 

Hi  is  the  horizontal  distance  in  inches  from  the  shoreline 
to  the  point  in  the  model  at  which  fresh  water  is  being 
extracted. 

Q2  is  the  rate  in  cubic  centimeters  per  minute  at  which 
fresh  water  is  being  extracted. 

Run  No.  1 

In  this  test  the  sea  chamber  of  the  small  model  was 
kept  full  of  salt  water  of  specific  gravity  1.040,  and  the 
salt  water  allowed  to  percolate  into  the  sand  without  any 
addition  of  fresh  water,  until  the  level  of  salt  water  in 
the  sand  reached  the  level  of  water  in  the  sea  chamber. 

Fresh  water  at  the  rate  of  75  cubic  centimeters  per  minute 
was  applied  at  H  =  10  inches.  Figure  14  in  Appendix  B  shows 
the  positions  of  the  interface  at  the  times  shown. 

Run  No.  2  to  Run  No.  13 

In  each  of  these  runs  an  equilibrium  position  was 
first  established  for  a  flow  rate  of  75  cc/min  of  fresh 
water  into  the  model.  The  fresh  water  was  put  in  at  a 
point  21  inches  from  the  shoreline.  The  outflow  rate  Q2 
was  varied,  as  well  as  the  distance  H  ^  ,  which  is  the  distance 
from  the  shoreline  of  the  point  of  extraction  of  the  fresh 
water.  Figures  15-26  in  Appendix  B  give  the  values  of 
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Q, Q2, H  and  for  each  run,  and  show  the  positions  of  the 
interface.  The  time  marked  on  each  position  of  the  inter¬ 
face  relates  to  the  time  which  had  elapsed  since  the 
commencement  of  pumping.  Some  of  the  Figures  show  the 
traces  made  by  dye  which  had  been  injected  at  various  points 
in  the  sand  while  extraction  was  taking  place.  The  density 
of  the  salt  water  for  all  runs  was  1.040. 

R  un  No.  14 

The  equilibrium  position  of  the  interface  for  an 
inflow  rate  of  75  cc/min  of  fresh  water  into  the  model  was 
established,  the  point  of  application  of  the  fresh  water 
being  21  inches  from  the  shoreline.  Fresh  water  was  then 
extracted  simultaneously  at  equal  rates  from  three 
extraction  points  placed  5,  15  and  25  inches  respectively 

from  the  shoreline.  The  total  rate  of  extraction  was 
50  cc/min.  Figure  27  in  Appendix  B  shows  the  positions  of 
the  interface,  after  the  elapsed  times  since  the  commence¬ 
ment  of  pumping  as  shown. 

Run  No.  15 

This  run  differed  only  from  Run  14  in  that  the  total 
extraction  rate  was  75  cc/min. 

Run  No.  16 

The  run  was  very  similar  to  Runs  No.  2  to  No.  13. 

The  point  of  input  of  fresh  water  into  the  system  was  26 
inches,  Q  =  75  cc/min,  Q2  =  50  cc/min.  and  H ^  = 


15  inches . 


(■;  i 
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Runs  No.  17  to  No.  19 

With  the  equilibrium  positions*  of  the  interface  as 
shown  at  A  of  Figure  28  in  Appendix  B,  the  extraction 

of  fresh  water  from  the  mo  del  was  discontinued,  while  the 
inflow  of  fresh  water  at  the  rate  of  75  cc/min.  was  main¬ 
tained  with  the  point  of  inflow  21  inches  from  the  shore¬ 
line.  The  Figure  30  shows  the  rate  of  recovery  of  the 
system,  that  is  the  rate  of  displacement  downwards  of  the 
interface  by  the  inflowing  fresh  water. 


CHAPTER  VII 


DISCUSSION  ON  DATA  OBTAINED  FROM 
LABORATORY  EXPERIMENTS 

The  laboratory  experiments  yielded  valuable  infor¬ 
mation  on  the  limitations  of  the  application  of  the  Ghyben- 
Herzberg  relation  to  a  coastal  aquifer  and  these  results 
and  limitations  are  discussed  here. 

A  set  of  curves  have  been  derived  which  relate  the 
various  parameters  which  manifest  themselves  in  the 
hydrostatic  balance  of  the  two  density  system. 

7  .  1  Discussion  on  Quant  it  at ive  Sect  ion  of  Lab oratory 
Studies 

Figures  10,  11  in  Appendix  B  show  a  set  of  curves 

which  have  been  obtained  as  discussed  in  Section  6.4,  and 
which  show  a  relation  between  q,  H,  L,  X,  Y  and  K  for  values 
obtained  in  both  models.  Figure  10  is  for  the  case  where 
H/L  =  0.93  while  Figure  11  is  for  the  case  where  H/L  =  0.43, 
and  H/L  =  0.54  in  the  small  and  big  models  respectively. 

From  each  q/KH  curve  at  a  particular  vertical  section, 
values  have  been  read  off  on  the  Y/H  axis,  where  a  vertical 
line  intersects  the  line  representing  q/KH.  Tables  III,  IV 
show  the  values  of  q/KH  and  the  corresponding  values  of 
Y/H. 
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T  AB  LE  III 

RELATION  OF  q/KH  AND  Y/H 
TAKEN  FROM  FIGURE  10 


q/KH  x  103 

Y/H 

.  9 

.  115 

1 .  25 

.  136 

1.35 

.  145 

1 .63 

.  155 

2 . 04 

.  175 

2  .  74 

.  205 

3  .  12 

.215 

3  .  70 

.  235 

<r 

00 

<r 

.280 

6 .04 

.330 

TABLE  IV 

RELATION  OF  q/KH  AND  Y/H 
TAKEN  FROM  FIGURE  11 

q/KH  x  103  Y/H 


1 . 0 

.  125 

1.5 

.  15 

1.87 

.  163 

2.34 

.  185 

3  .  24 

.  215 

4 . 9 

.26  0 

5.62 


.295 
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Figures  12,  13  in  Appendix  B  show  the  relation  by  means  of 

a  plot  of  Y/H  versus  q/KH.  The  lines  A  t o * H  rn  Figure 

10  plotted  on  double  logarithmic  paper  have  a  slope  of 
one-half,  and  hence  indicate  the  parabolic  shape  of  the  inter¬ 
face.  Lines  J  and  K  both  possess  a  slope  less  than  one- 

half  ,  and  when  the  q/KH  values  associated  with  these  two 
curves  are  considered,  it  seems  to  indicate  that  there  is 
a  limiting  value  of  q/KH,  in  the  order  of  4  x  10-^,  beyond 
which  the  parabolic  relation  of  the  shape  of  the  interface, 
does  not  hold  for  this  particular  H/L  ratio.  In  Figure  11 
for  H/L  ratio  of  0.43  in  the  small  model  and  H/L  =  0.54 
in  the  big  model,  this  limiting  value  of  q/KH  has  been 
exceeded,  but  the  parabolic  relation  of  the  shape  of  the 
interface  still  holds.  A  comparison  of  Figures  10  and  11 
clearly  indicate  the  effect  of  the  distance  from  the  shore¬ 
line  of  the  point  of  concentration  of  inflow  of  the  fresh 
water  into  the  model.  In  Figure  11  a  distinct  break  from 
the  parabolic  relation  of  the  shape  of  the  interface  is 
clearly  indicated,  and  the  break  appears  to  centre  around 
an  X/H  value  of  1.0,  which  corresponds  to  the  point  of 
application  of  the  fresh  water  in  either  model.  Figure  10, 

11  give  some  indication  of  the  reliability  of  the  method 
of  treatment  of  the  data  in  the  attempt  to  secure  some 
dimensionless  curves  which  can  be  used  to  relate  the  various 
parameters  discussed.  The  results  of  the  big  model  are 
very  much  in  agreement  with  the  results  of  the  small,  model. 
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7 . 2  L imit  at  ions  o  f  Use  o  £  Curves  o  f  F igur  e  s  10 ,  1 1 

The  q/KH  and  Y/H  values  obtained  in  Figures  10,  11 

have  all  been  derived  from  model  tests  in  which  the  density 
of  salt  water  was  1.040.  To  make  use  of  the  curves  in 
Figures  12,  13,  the  assumption  has  to  be  made  that  Y  varies 


1 inear ly  with 


- ,  and  hence  the  values  of  Y  obtained  in 

P-R 


the  prototype  have  to  be  multiplied  by  a  factor  of  about 


25 

"411 


to  correspond  to  the  Y  of  the  curves.  With  this  assumption 
made,  and  provided  that  the  H/L  ratio  of  the  prototype 
corresponds  to  one  of  the  H/L  ratios  of  the  curves,  the  value 
of  q/KH  corresponding  to  the  corrected  value  of  Y/H  from  the 
prototype  can  be  obtained  from  either  Figure  12  or  Figure  13. 
When  this  value  has  been  obtained  there  are  then  only  two 
unknowns--q  and  K  of  the  prototype.  The  value  of  K  can  be 
determined  by  one  of  several  methods  including  the  use  of 
tracers  for  velocity  determinations.  The  usefulness  of  the 
curves  in  Figures  10,  11  is  further  enhanced,  by  the 

observation  that  if  the  Figures  10,  11  are  superimposed, 

there  is  a  close  fit  of  the  curves  for  the  portions  of  the 
curves  corresponding  to  X/H  less  than  1.0.  This  seems  to 
indicate  that  either  curve  can  be  used  for  H/L  values  ranging 
from  1.0  to  0.43  provided  the  portions  of  the  curve  for  which 
X/H  does  not  exceed  1.0  are  used. 

The  range  of  values  of  q/KH  for  which  curves  have 


been  obtained  is  very  small,  and  shows  clearly  the  rather 
stringent  limits  of  application  of  the  data.  With  the 
models  used,  it  was  impossible  to  secure  values  outside  of 
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the  range.  The  particular  prototype,  for  example,  in  which 
the  author  is  interested,  would  necessitate  some  very  small 
and  impractical  flow  rates  of  fresh  water  into  the  model, 
or  conversely  a  large  increase  in  the  values  of  K  and  H 
used.  Large  increases  in  the  value  of  the  permeability  of 
the  medium  used  in  the  model  will  certainly  present 
difficulties  because  of  the  possibility  of  a  change  over 
from  laminar  to  turbulent  flow.  The  other  alternative  then, 
is  a  sufficiently  long  model  so  that  the  value  of  H  can  be 
substantially  increased. 

7 . 3  P  er meab i 1  it  y  Values  Ob  t  a  in  ed  in  Different  Tests 

The  values  of  permeability  found  in  the  models  for 
different  flow  rates  are  as  follows: 

Big  Mo  del .  H  =  114  inches 


Q 

in  c  c  /  m  in  . 

K 

x  10  ^  cms/sec. 

7  .  3 

2.03 

12 . 5 

2 . 5 

14 

6.02 

62 

6.33 

B  ig 

Model.  H  = 

66 

inches 

Q 

in  c  c  /  m  in  . 

K 

x  10  ^  cms/sec. 

14 

6 .03 

30 

8.07 

58 

11.7 

80 

9.3 

t 


Big  Model 


Permeability  from  Permeameter  test  K.  =  2,68  x  10  ^ 
cms/sec. 

From  velocity  measurement  for  Q  =  14  c.  c/m  in. 

K  =  7.02  x  10“ “  cms/sec. 

Small  Model  for  H  =  26  inches 


Q 

in  cc/min . 

K 

cms/sec. 

12 

.  2  7 

40 

.  24 

80 

.36 

Sma  1 1 

Mo  del  for 

H  = 

12  inches 

Q 

in  cc/min. 

K 

cm/ s  ec . 

5 

.  30 

10 

.  33 

20 

.  25 

40 

.32 

Sma 1 1  Mode  1 

From  Permeameter  test  K  =  1.68  x  10 “2  cms/sec. 

From  velocity  measurement  Q  =  16  cc/min.  K  =  .875  cm. /sec. 

The  results  obtained  for  tests  in  the  small  model  tend  to 
show  a  greater  degree  of  consistency  for  the  values  of 
permeability  for  different  flow  rates  than  in  the  big  model. 

Yet  in  the  case  of  the  big  model,  there  is  a  closer  fit  of 
values  of  permeability  obtained  in  the  model  and  from  the 
permeameter  and  velocity  methods.  The  smaller  model  on  the 
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other  hand  yielded  values  of  permeability  which  were  smaller 
than  those  obtained  by  a  velocity  measurement,  and  are 
almost  twenty  times  bigger  than  the  value  of  permeability 
obtained  in  the  permeameter  test.  It  can  only  be  suggested 
as  an  explanation  of  this,  that  the  packing  employed  and 
hence  the  void  ratios  in  the  two  tests,  that  is,  in  the  small 
model  and  in  the  permeameter  were  entirely  different. 

Tables  I,  II  in  section  6.6  show  the  velocities  of  flow 
obtained  by  injecting  dye  at  varying  depths  which  may  explain 
the  consistency  of  the  values  of  permeability  obtained  in 
the  small  model  and  the  variation  with  flow  rate  and  hence 
with  depth  in  the  model,  of  values  obtained  in  the  big  model. 

7 . 4  General  Obs  er vat  ion  on  Quant  it  at ive  Sect  ion  o  f  Research 

The  observation  was  made  that  in  one  or  two  cases, 
there  appeared  to  be  a  discrepancy  in  the  positions  of  the 
interface  for  certain  flow  rates.  For  example,  flow  rates  of 
7.3  and  12.5  cubic  cms .  per  minute  yielded  values  of 
permeability  which  are  close  to  each  other,  but  which  in  turn 
are  about  one  third  the  values  obtained  in  the  other  tests 
under  similar  conditions  of  the  point  of  inflow. 

An  attempt  was  made  to  reproduce  the  flow  of  12.5  cub  ip 
centimeters  per  minute,  but  due  to  a  slight  error  of  judgement 
a  flow  of  14  cc/min.  was  established.  Tables  XVIII,  XXXII 
show  a  change  in  the  position  of  the  interface  in  the  two  tests, 
in  spite  of  the  small  difference  in  the  flow  rates.  Yet  in 
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the  final  analysis  the  q/KH  values  obtained  from  the  tests 
of  the  two  flow  rates  are  in  good  agreement,  and  fit  very 
well  the  curves  of  Figure  12. 

This  seems  to  indicate  that  the  permeability  of  the 
model,  in  this  particular  case  the  big  model,  must  have 
changed  somewhat  during  the  period  between  the  two  tests. 
This  is  further  indicated  by  two  positions  of  the  inter¬ 
face  obtained  for  the  same  flow  rates  of  33  centimeters  per 
minute  and  are  shown  in  Figure  6.  Whereas  a  small  change  in 
sectional  width  of  the  model  due  to  bulging  by  virtue  of 
the  weight  of  the  sand  might  have  affected  the  overall 
values  of  permeabilities  discussed  in  section  6.5,  it  is 
doubtful  that  this  change  in  section  alone  can  have  such 
large  effects.  The  question  then  arises  as  to  the 
possibility  of  changes  in  the  degree  of  saturation  of  the 
sand,  with  possible  formation  of  air  pockets  in  the  models, 
during  the  change-over  from  one  test  to  the  other.  The 
results  definitely  cast  doubts  on  the  reliability  of 
accepting  a  certain  position  of  the  interface  as  indicative 
of  the  rate  of  inflow  of  fresh  water  into  the  aquifer,  and 
that  one  should  rather  e mp loy  the  whole  process  of  the 
treatment  of  the  data  as  derived  in  this  work,  and  make 
use  of  the  q/KH.  versus  Y/H  relation. 

7  .  5  The  Gh  yben -Her  z  berg  Relation  in  _a  Non-P  ump  ing 
System 


The  second  and  fourth  columns  of  Tables  V  to  XXII  in 
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Appendix  A  show  the  values  of  the  depths  from  mean  sea  level 
to  the  interface  as  actually  measured  in  the  models,  as 
compared  with  values  calculated  from  an  application  of  the 
Ghyb en-Her zb er g  relation.  Tables  V  to  XIV  relate  to  the 
small  model  while  Tables  XV  to  XXII  are  for  the  big  model. 

In  each  case  the  flow  rates  of  fresh  water  into  the  model, 
and  the  distance  from  the  shoreline  of  the  point  of 
application  of  the  fresh  water  are  given. 

In  the  case  of  the  small  model,  it  was  impossible  to 
measure  the  height  of  the  water  table  above  mean  sea 
level  ( Y 1 )  to  an  accuracy  greater  than  one-sixteenth  of  an 
inch.  Since  this  value  of  Y  ^  has  to  be  multiplied  by 


P„ 


P-fi. 


=  25  for  the  theoretical  depth  to  the  inter¬ 


face  according  to  the  Ghyben-Her zberg  relation,  then  a 
variation  of  one  to  two  inches  from  the  theoretical  depth 
should  not  be  considered  to  be  significant.  It  was  possible 
to  secure  a  greater  degree  of  accuracy  of  measurement  in 
the  big  model,  which  might  account  for  the  closer  fit  of 
the  values  under  consideration. 

Consideration  of  the  Tables  V  to  XXII  show” 

(a)  Very  near  to  the  shoreline  the  Ghyb en-Her zb er g 

relation  does  not  hold.  In  some  cases,  the  calculated 
values  are  less  than,  and  in  other  cases  greater 
than,  the  measured  values  of  the  depth  to  the  inter¬ 
face.  The  discrepancy  in  the  values  of  the  depth  to 
the  interface  was  expected  however,,  since  several  of 
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the  Figures  show  a  large  seepage  area  at  the  shore¬ 
line  through  which  fresh  water  is  flowing.  Other 
Figures  show  a  concentration  of  curved  flow  lines  in 
this  area  with  subsequent  increases  in  velocities  of 
flow  and  turbulent  flow  in  the  region.  The  deri¬ 
vation  of  the  Ghyben-Her zberg  relation  assumed  a 
statical  system, and  hence  is  not  applicable  in  this 
region  near  the  shoreline. 

(b)  The  distance  (H)  from  the  shoreline  of  the  point  of 
application  of  the  fresh  water  must  be  taken  into 
account  in  a  consideration  of  the  reliability  of  the 
Ghyben-Her zberg  relation.  This  is  indicated  in  the 
case  of  the  small  model ,  but  is  more  clearly  seen  by 
a  consideration  of  Tables  V  to  XI  and  Tables  XIX  to 
XXII  for  the  big  model.  In  Tables  XV  to  XVIII  when 
the  fresh  water  was  applied  at  a  point  114  inches 
from  the  shoreline,  the  measured  values  and  the 
calculated  values  of  the  depth  to  the  interface  are 
not  significantly  different.  But  in  the  tests  in 
which  the  distance  from  the  shoreline  of  the  point  of 
application  of  the  fresh  water  is  66  inches,  then 
there  are  differences  in  the  values  of  the  depth  to 
the  interface  as  compared  to  the  calculated  values 
(Tables  XIX  to  XXII).  While  there  are  large 
differences  in  the  values  seaward  of  the  point  of 
application  of  the  fresh  water,  the  differences  are 
even  larger  landwards  of  this  point,  and  in  neither 


case  does  the  calculated  values  from  a  consideration 


of  the  Ghyben-Her zb  erg  relation  approach  the  measured 
values.  There  appears  to  be  a  ponding  up  of  the 
water  in  the  saturated  zone  above  mean  sea  level  land¬ 
wards  of  the  point  of  application  of  the  fresh  water, 
which  then  yields  calculated  values  of  the  depth  to 
the  interface  which  in  some  cases  are  fifty  per  cent 
too  large.  (Tables  XIX  to  XXII). 

This  then  means,  that  as  the  H/L  ratio  decreases, 
so  does  the  reliance  which  can  be  placed  on  the  values 
obtained  for  the  depth  to  the  interface,  from  a 
consideration  of  the  Ghyben-Her zberg  relation. 

7 . 6  Cons  id  e  rat  ion  o  f  the  Ghyb en-Her zberg  R  e 1 at  ion  in  a 
P  umped  System 

Figures  20-34  show  the  various  positions  of  the  inter¬ 
face  with  time  under  conditions  of  extraction  of  fresh  water 
from  various  point  sources  in  the  small  model.  Section  6.8 
gives  a  description  of  the  runs  corresponding  to  Figures  20- 
34.  Some  of  the  Figures  also  show  the  position  of  the  water 
table  before  pumping  commenced,  and  the  position  of  the 
water  table  just  before  the  cessation  of  pumping  when  no 
further  change  in  the  position  of  the  interface  occurred. 

These  two  curves  are  labelled  A  and  R  respectively  in  the 
upper  diagrams  of  Figures  15  to  29  in  Appendix  B.  Position 
B  of  the  water  table  was  obtained  in  all  cases  within  five 
minutes  of  the  commencement  of  pumping.  Inspection  of 
Figures  15  to  29  indicate  that  except  in  the  cases  where  a 
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cone  of  salt  water  is  drawn  into  the  well,  and  within  the 
previously  mentioned  limits  of  the  accuracy  of  measurement 
in  the  small  model,  a  reasonably  close  value  is  obtained  to 
the  measured  value  of  the  initial  and  final  position  of  the 
interface  under  pumping  conditions,  when  the  Ghyben- 
Herzberg  relation  is  applied  to  curves  A  and  B,  the  initial 
and  final  positions  of  the  water  table.  Since,  however, 
the  ratio  of  the  times  for  the  position  of  the  water  table 
and  the  interface  to  reach  equilibrium  positions  is  of  the 
order  of  five  minutes  to  five  hours,  or  1:60,  it  would 
appear  that  application  of  the  Ghyb en-Her zb er g  relation  to 
find  the  position  of  the  interface  during  the  periods 
between  the  two  equilibrium  positions  is  indeed  very  subject 
to  error. 

Superimposed  on  the  effect  of  the  distance  from  the 
shoreline  of  the  point  of  application  of  the  fresh  water 
as  discussed  in  section  7.5,  is  the  effect  of  the  distance 
of  the  pumping  source  from  the  shoreline.  Figure  15  Run 
No.  2  is  a  good  indication  of  this.  The  combination  of  a. 
pumping  source  nineteen  inches  from  the  shoreline,  with  a 
point  of  concentration  of  inflow  of  fresh  water  twenty-one 
inches  from  the  shoreline,  had  the  effect  of  producing  a 
change  in  the  water  table  only  in  the  region  landwards  of  a 
point  fifteen  inches  from  the  shoreline  in  the  model,  yet 
with  the  whole  interface  shifting  its  position  as  indicated 


in  the  F igur  e  15 . 
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Figures  19,  20,  27  Runs  6,  7  and  14,  show  the  effect 

of  ponding  up  of  water  in  the  saturated  zone  above  sea  level 
near  the  point  of  inflow  of  fresh  water  into  the  aquifer, 
and  the  Ghyben-Her zberg  relation  cannot  be  used. 

In  the  case  where  a  cone  of  salt  water  extends 
up  into  the  well,  the  Ghyben-Her zberg  relation  is  no  longer 
applicable.  This  is  indicated  by  Figure  2,0  Run  No.  7  in  which, 
were  the  Ghyben-Her zberg  relation  applied  to  the  point  nine 
and  one-half  inches  from  the  shoreline, then  the  apex  of  the 
cone  should  have  been  located  at  a  depth  of  four  to  five 
inches  instead  of  the  measured  value  of  one  and  one-half 
inches.  This  breaking  down  of  the  relation  is  even  more 
marked  in  the  case  of  Run  No.  9  Figure  22,  where  the 
position  of  the  apex  of  the  cone  differs  from  the  position 
it  should  occupy  according  to  the  Ghyben-Her zberg  relation 
by  as  much  as  four  inches  in  the  model.  Run  No.  11  Figure 
24  further  indicates  the  unreliability  of  the  Ghyben- 
Herzberg  relation  when  applied  to  a  case  where  salt  water 
has  intruded  into  the  pumping  source. 

7  7  Gener  a  1  Ob  s  er  vat  ions  From  the.  P  u  mp  e  d  Tests 

Within  the  limits  of  a  two  dimensional  model  the 
following  observations  are  made: 

(a)  Except  for  a  region  very  close  to  the  shoreline,  and 
with  the  bottom  of  the  well  at  or  near  sea  level, 
there  appears  to  be  a  safe  extraction  rate  of  fresh 
water  from  the  aquifer,  which  can  be  maintained 


without  the  danger  of  drawing  salt  water  into  the  well, 
irrespective  of  the  distance  of  the  well  from  the 
shoreline.  This  extraction  rate  appears  to  be  of  the 
order  of  seventy  per  cent  of  the  total  inflow  of 
fresh  water  into  a  region  of  the  aquifer. 

(b)  In  Run  No.  15  Figure  28,  an  extraction  rate  of 

25  cc/min,  was  maintained  simultaneously  in  each  of 
three  wells  placed  at  5,  15  and  25  inches  respectively 

from  the  shoreline.  Figure  28  shows  the  different 
locations  of  the  interface  at  various  elapsed  times 
since  the  commencement  of  pumping.  The  Figure  shows 
that  after  three  hours,  the  well  placed  five  inches 
from  the  shoreline  was  drawing  salt  water,  while  the 
other  wells  were  unaffected  by  the  salt  water,  and 
that  after  five  and  one-half  hours,  when  the  inter¬ 
face  had  reached  an  equilibrium  position,  the  wells, 
placed  at  five  and  fifteen  inches  respectively  from  the 
shoreline  were  drawing  salt  water,  while  the  well  at 
twenty-five  inches  was  unaffected.  The  test,  then, 
demonstrates  that  the  wells  closer  to  the  shoreline, 
tend  to  become  contaminated  with  sea  water  before  the 
wells  which  are  further  inland,  and  that  even  by 
splitting  the  flow  rates  between  a  number  of  wells, 
which  all  lie  on  the  shortest  line  between  the  point 
of  the  inflow  of  the  fresh  water  and  the  shore  line. 


the  total  extraction  rate  of  fresh  water  must  be  less 


than  the  flow  rate  of  fresh  water  into  the  aquifer, 
without  the  danger  of  salt  water  intrusion  into  the 


wells . 

(c)  Figure  3  in  section  3.2  shows  the  cone  of  salt  water 

extending  upwards  with  its  axis  vertical  towards  the 
point  of  extraction.  This  is  usually  the  way  the  cone 
of  salt  water  is  shown  in  textbooks,  (Ref  10)  and  can 
easily  be  interpreted  to  be  an  indication  that  the 
axis  of  the  cone  is  always  vertical.  Figures  20,  22 

show  that  this  is  not  necessarily  the  case,  but  that 
the  axis  of  the  cone  of  indrawn  sea  water  can  lie 

and  generally  does  lie  at  an  angle  inclined  to  the 
vertical,  due  to  the  prevalent  direction  of  flow  of 
fresh  water  as  it  is  drawn  towards  the  well.  Figure 
26  Run  No.  13  shows  the  cone  of  salt  water  being 
formed  on  the  landward  side  of  the  point  of 
extraction  of  the  fresh  water.  The  importance  of 
these  phenomena  is  easily  recognized  when  one  is 
considering  the  placement  of  boreholes  in  the  field  to 
check  on  salt  water  intrusion  into  the  well. 

( d)  The  paths  traced  by  dye  which  had  been  injected  at 
various  points  in  the  model  are  shown  in  several  of 
the  Figures  15  to  22.  The  traces  show  that  water  in 
all  horizontal  sections  of  the  aquifer  are  in  motion, 
with  the  velocities  of  flow  having  a  larger  vertical 
downward  component  near  the  points  at  which  water  is 
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put  into  the  model.  In  horizontal  sections  near  the 
sea  level  line,  the  flow  is  almost  horizontal,  and 
for  deeper  sections  the  flow  paths  follow  more  closely 
the  contour  of  the  interface.  The  dye  traces  clearly 
show  the  flow  of  water  which  bypasses  the  well,  and 
also  indicate  the  loss  of  this  water  to  the  sea, 
through  the  seepage  face  near  the  shoreline. 


CHAPTER  VIII 


APPLICATION  OF  RESULTS  OF  LABORATORY  STUDIES 
TO  GROUND  WATER  PROBLEM  OF  BARBADOS 

The  chapter  serves  as  a  basis  of  discussing  the 
usefulness  of  the  research  and  an  attempt  is  made  very 
briefly  to  apply  some  of  the  results  and  observations  from 
the  laboratory  work  to  the  conditions  in  Barbados,, 

8,1  Previous  Ideas  on  Mode  of  Occurrence  of  Ground  Water 
in  B  arb  ados 

Figure  35  in  Appendix  B  taken  from  a  report  by  S enn 
(Ref  16)  shows  the  views  previously  held  by  that  author  on 
the  mode  of  occurrence  of  ground  water  in  the  coastal  aquifer, 
S enn  shows  fresh  water  occupying  the  entire  volume  of  the 
region  bounded  by  the  sea  at  the  shoreline;,  the.  clay-coral 
interface  at  the  bottom,  and  on  the  landward  side  by  the 
sloping  face  of  the  clay  interface.  This  could  have  been 
correct,  had  the  depth  of  coral  within  the  aquifer  been  of 
a  much  lesser  value  and  would  have  prohibited  the  stabili¬ 
zation  described  by  the  Ghyben-Her zberg  relation.  But  with 
the  thickness  of  coral  exceeding  two  hundred  feet  in  this 
region,  and  with  the  height  above  sea  level  of  the  water 
table  not  exceeding  three  feet  in  the  coastal  aquifer,  the 
ground  water  must  exist  to  some  degree  in  the  form  prescribed 
by  the.  Ghyben-Her  zberg  relation,  that  is,  that  the  fresh 


water  contained  in  the  aquifer  must  everywhere  in  the  aquifer 
be  underlain  by  sea  water,  except  for  regions  along  the  land¬ 
ward  boundary  of  the  sheetwater  where  the  coral  is  of  in¬ 
sufficient  thickness.  Subsequently  the  volume  in  storage  as 
shown  in  Figure  35  is  incorrect  and  should  be  very  nearly 
halved  for  a  truer  value.  Figure  39  in  Appendix  B,  referred 
to  in  section  4.2  gives  the  correct  picture  of  the  mode  of 
occurrence  of  the  ground  water  in  this  coastal  aquifer. 

8 . 2  P  r  ob lems  Associated  With  the  Mis  interpret  at  ion  o  f  the 

0 c.  cur  r  enc  e  o  f  the  System 

As  discussed  in  section  8.1,  the  opinion  existed  that 
the  ground  water  in  the  coastal  aquifer,  shown  as  sheetwater 
in  Figure  33  in  Appendix  B,  existed  as  a  rather  static  lake 
and  that  extraction  from  a  particular  source  drew  from  the 
lake  without  any  danger  of  salt  water  intrusion.  This  has 
been  proven  to  be  incorrect  as  discussed  in  section  8.1. 

There  is  the  danger  of  salt,  water  intrusion  and  considerable 
loss  of  fresh  water  to  the  sea  along  the  seepage  faces  near 
the  shoreline  as  has  been  demonstrated  in  the  models. 

8 . 3  Suggested  R  erne  dial  Method  s 

(a)  The  practice  of  withdrawing  vast  quantities  of  fresh 
water  from  a  single  source  near  the  shoreline  has 
been  shown  to  present  considerable  danger  of  salt 
water  intrusion,  and  hence  should  be  avoided.  It 
will  be  a  far  better  practice  to  extract  smaller 
quantities  from  several  points  within  the  coastal 
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aquifer,  but  as  far  inland  as  is  practicable.  This 
practice  will  reduce  the  danger  of  salt  water  con¬ 
tamination,  and  at  the  same  time  prevent  much  of  the 
loss  of  the  fresh  water  to  the  sea  which  has  been 
shown  in  the  models  to  be  continually  taking  place 
because  of  the  dynamic,  nature  of  the  system. 

(b)  Due  to  the  unavailability  of  reliable  data  as  to  the 

permeability  of  the  coral,  no  attempt  can  be  made 
here  to  make  use  of  the  quantitative  aspects  of 
Figures  10,  11.  When  these  data  are  available,  and 

provided  that  the  flow  rates  are  within  the  range  of 
the  model  studies,  values  of  the  flow  rate  of  fresh 
water  per  foot  width,  of  aquifer  can  be  calculated 
and  will  in  deed  be  very  useful.  If  the  range  of  the 
model  studies  of  the  present  work  prove  to  be  too 
limited,  then  the  suggestion  discussed  in  Section  7.2 
can  be  implemented  and  a  longer  model  built  and 
studied. 

(c)  Because  of  the  fallacy  that  the  ground  water  exists 
as  a  common  static  lake  discussed  in  Section  8.2,  the 
ground  water  along  the  south  coast  of  the  island  has 
not  been  developed  and  should  certainly  be  developed. 
The  flow  of  fresh  water  to  the  sea  in  the  area  is 
certainly  taking  place  and  this  fresh  water  is  being 
lost. 

(d)  Whereas  the  application  of  the  Ghyb en-Her zb er g  relation 
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can  at  least  give  some  ideas  on  the  volume  of  fresh 
water  in  storage  by  a  knowledge  of  the  approximate 
position  of  the  interface  in  a  non-pump ed  system, 
not  much  reliability  should  be  placed  on  values 
obtained  by  the  application  of  the  Ghyb en-Her zb er g 
relation  on  a  pumped  system.  For  more  reliable- 
results,  it  should  first  be  ascertained  that  the 
system  has  reached  an  equilibrium,  and  that  the  lag 
effect  as  discussed  in  Section  7.6  is  not  in 
existence.  It  also  follows  from  the  discussion  in 
Section  7.7  that  a  probe  in  one  borehole  is  in¬ 
sufficient  to  establish  the  absence  of  a  cone  of 
indrawn  salt  water. 

There  are  obvious  other  suggestions  which  could  be 
made,  but  which  do  not  follow  as  a  direct  consequence  of  this 


study , 


CHAPTER  IX 


SUMMARY  AN B  CONCLUSIONS 

The  chapter  summarizes  the  findings  of  the  laboratory 
research,  and  in  it  also  the  author  makes  some  suggestions 
to  further  research  of  the  topic. 

9 .  1  S  umma r y 

Previous  literature  on  the  subject  shows  that  only 
general  qualitative  studies  of  the.  Gh  yb  en-H  er  zb  er  g  relation 
have  been  investigated.  No  attempt  has  previously  been  made 
to  secure  a  set  of  universal  curves,  which  relate  the  various 
parameters  associated  with  such  a  two  density  system.  The 
present  study,  then,  provides  such  a  set  of  curves,  although 
admittedly  they  are  rather  limited  in  their  scope  of 
application.  The  use  of  two  models,  one  serving  as  a 
prototype  of  the  other,  provided  a  means  of  checking  the 
results,  and  proves  the  feasibility  of  such  quantitative 
studies  of  the  two  density  system  in  a  coastal  aquifer. 

The  limitations  on  the  use  of  the  Ghyb en-Her zber g 
relation  have  been  observed  and  discussed,  and  are  especially 
marked  in  the  case  of  a  pumped  aquifer.  Certain  quantitative 
data  have  been  collected  relating  the  positions  of  the  inter¬ 
face  with  elapsed  periods  of  time  after  the  commencement  of 
pumping.  It  is  hoped  that  with  additional  data  from  an 
actual  prototype,  some  relation  between  the  position  of  the 
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interface,  rate  of  pumping  and  position  of  the  interface  can 
be  established. 

9 . 2  C  one  lus ions 

A,  Application  of  the  Ghyben-Her zb er g  relation  to  a  non- 
pumped  aquifer,  or  in  cases  where  the  rate  of  pumping 
is  small  enough  that  the  system  can  still  be  con¬ 
sidered  a  statical  one,  yields  values  of  the  depth 

to  the  interface  which  for  practical  purposes  are 
acceptable,  except  for  cases  mentioned  in  B. 

B,  The  distance  from  the  shoreline  of  the  point  of 
inflow  of  fresh  water  into  the  aquifer,  relative  to 
the  overall  length  of  the  per me able  stratum  in  the 
aquifer  is  very  important,  and  has  considerable 
effect  on  the  application  of  the  Ghybea-Herzberg 
relation.  It  is  safe  to  apply  the  Ghyben-Her zberg 
relation  only  in  the  region  seawards  of  the  point  of 
inflow  of  fresh  water  into  the  aquifer.  As  the  ratio 
of  the  distance  from  the  shoreline  of  the  inflow  of 
fresh  water  to  the  length  of  the  permeable  stratum 
decreases,  so  does  the  reliability  which  can  be 
placed  on  the  results  obtained  by  an  application  of 
the  Ghyben-Her zberg  relation, 

C ,  In  systems  where  salt  water  has  been  drawn  into  the 
well,  the  Ghyben-Her zberg  relation  cannot  be  used. 

In  most  cases,  the  Ghyben-Her zberg  relation  gives  a 
depth  to  the  interface,  which  can  be  three  or  four 
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t ime  s  too  large, 

D.  The  dynamical  nature  of  the  two  density  system  has 
been  established,  A  hydraulic  gradient  of  the  water 
table  roust  be  maintained,  with  a  subsequent  continuous 
f low  o  f  fresh  water  t  o  the  sea, 

E.  There  is  not  necessarily  a  unique  position  of  the 
interface  for  a  particular  flow  rate  of  fresh  water 
into  the  aquifer.  The  q/KH  value  for  a  particular 
flow  rate  into  the  aquifer  is  unique  however, and  is 
the  criterion  that  should  be  established  and  used, 

F,  Figures  10,  11  provide  a  universal  set  of  curves  which 
can  be  used  to  establish  relations  of  the  depth  to 

the  interface,  permeability  of  the  model,  and  the 
distance  from  the  shoreline  of  the  point  of  application 
of  the  fresh  water, 

G,  There  appears  to  be  a.  range  of  values  of  the  rate  of 
inflow  of  fresh  water  into  the  aquifer,  beyond  which 
the  parabolic  relation  of  the  shape  of  the  interface 
does  not  hold.  It  is  believed  that  this  range  is 
dictated  by  the  change-over  from  laminar  to  turbulent 
flow  between  the  soil  grains  of  the  permeable  medium, 

9 , 3  Recommendat ions  for  Future  Study 

A,  For  the  particular  prototype  in  which  the  author  was 
interested,  a  two  dimensional  model  was  considered 
to  be  adequate  in  order  to  obtain  a  first  approach 
to  the  problem.  A  research  program  involving  the  use 


of  a  three  dimensional  model  seems  then  to  suggest 
itself.  It  would  also  be  desirable  to  make  use  of 
the  three  dimensional  electric  analogue  to  verify 
streamline  patterns. 

Studies  should  be  undertaken  of  cases  where  layers  of 
impermeable  material  are  placed  in  the  permeable 
medium,  in  vertical  and  horizontal  layers  respectively. 
The  effect  of  tides  on  such  a  two  density  system,  has 
not  been  taken  into  consideration.  The  model  can 
easily  be  adapted  to  such  a  study,  and  it  is  suggested 
that  this  effect  be  studied  in  a  future  research 
program. 

It  has  been  suggested  that  one  of  the  causes  of  the 
change  in  permeability  of  the  mode  1 ,  as  experienced 
in  the  various  runs,  could  have  been  due  to  some 
chemical  action  between  the  salt  water  and  the  sand. 
This  chemical  action  could  be  in  the  form  of  a 
cementing  material  being  deposited,  with  a  subsequent 
change  in  the  size  of  the  porespace.  The  suggestion 
should  be  investigated. 


B  IB  BIOGRAPHY 


66 


1.  Badon  Ghyben,  W.,  Nota  in  verband  met  de  voor genomen 

putbor ing  nabij  Armsterdam  (Notes  on  the  probable 
results  of  the  proposed  well  drilling  near 
Amsterdam) .  Tijdschrift  van  het  Koninklijk 
Institut  van  Ingenieurs,  The  Hague,  p.  21,  1888- 

1889  . 

2.  California  Department  of  Public  Works,  Progress  report 

on  investigation  of  use  of  earthen  cut-off  walls 
for  prevention  of  sea  water  intrusion;  Division 
of  Water  Resources,  Memo.,  February  29,  1952. 

3.  Ohrt,  F.,  Water  development  and  salt  water  intrusion 

on  Pacific  Islands;  Jour , Amer .  Water  Works  Assoc., 
v  o 1 .  39,  pp.  979-988  ,  1947. 

4.  Riddel,  J.O.,  Excluding  salt  water  from  island  wells-- 

a  theory  of  the  occurrence  of  ground  wat  er  based 
on  experience  at  Nassau,  Bahama  Islands,  Civil 
Eng.,  vol.  3,  pp .  383-385,  1933. 

5.  Martin-Kaye,  P . ,  Extracts  from  Progress  Report  No. 11 

Geological  Surveys  Windward  Islands,  St.  Lucia, 
February,  1961. 

6.  Herzberg,  B»,  Die  Was s er ver s or gung  einiger  Nordeubader, 

(The  water  supply  on  parts  of  the  North  Sea 
coast);  Jour,  gasb e leucht ung  and  Wa s s er v er s or  gun g 
vol. 44,  pp.  815-819,  Munich,  1901. 

7.  Pennink,  Investigations  for  ground  water 

supplies;  Trans.  Amer.  Soc,  Civil  Engrs . ,  vol.  54-D, 
pp.  169-181,  1905. 

8.  d'Antrimont,  R.,  Preliminary  notes  upon  a  new  method 

of  studying  experimentally  the  phenomena  of  ground 
water  in  freely  pervious  aquifers,  Soc.  Geo  1 . 
Belgique  Ann ales,  vol. 32,  pp.  M  115-  M  120,  1905. 

9.  Nomitsu,  T.,  Y.  Toyohara,  and  R.  Kamimoto,  On  the  contact 

surface  of  fresh-  and  salt-water  near  a  sandy  sea¬ 
shore,  Mem. College  Sci. ,  Kyoto  Imp.  Univ., 

Ser.A,  vol. 10,  no. 7,  pp.  2  7  9-  302,  1927. 

10.  Todd,  D.K.,  Ground  water  Hydrology,  John  Wiley  and 
Sons,  New  York,  1960. 


¥  v 


* 


67 


11.  Banks,  H.O.,  and  R.C.  Richter,  Sea-water  intrusion  into 

ground  water  basins  bordering  the  California  coast 
and  inland  bays,  Trans.  Amer .  Geophysical  Union, 
v  o 1 .  34,  pp.  575-582,  1953. 

12.  Simpson,  T.R.,  Salinas  Basin  investigation,  Bull. 52,  Calif. 

Div.  Water  Resources,  Sacramento,  230  pp . ,  1946. 

13.  Banks,  H.O.,  and  Bookman,  M.,  Proposed  investigational 

work  for  control  and  prevention  of  sea  water 
intrusion  into  ground  water  basins;  Report  to  State 
Water  Resources  Board,  Div.  of  Water  Resources, 

State  of  California  33  pp . ,  Sacramento,  1951. 

14.  Laverty,  F.B.,  Jordan,  L.W.,  and  van  der  Goot,  H.A.; 

Report  on  tests  for  the  creation  of  fresh  water 
barriers  to  prevent  salinity  intrusion,  performed 
in  West  Coastal  Basin,  Los  Angeles  County, 

California,  70  pp . ,  1951. 

15.  Jacobs,  J.A.,  Russell,  R.D.,  and  Wilson,  J.T.,  Physics 

and  geology;  McGraw-Hill  Book  Company,  Inc. , 

Toronto,  p.  301,  1959. 

16.  Senn,  A.,  Report  of  the  British  Union  Oil  Company  Limited 

on  geological  investigations  of  the  ground  water 
resources  of  Barbados,  B.W.I.,  Barbados,  March,  1946. 

17.  Mart  in-Kaye ,  P.  ,  Barbados  Water  Supply;  Report  to  the 

Ministry  of  Natural  Resources  and  Agriculture,  The 
West  Indies,  July  1961. 

18.  Lamb,  T.W.,  Soil  Testing  for  Engineers;  John  Wiley  &  Sons, 

Inc. ,  New  York  1951. 

19.  Todd,  D.K.,  An  abstract  of  literature  pertaining  to  sea 

water  intrusion  and  its  control;  Tech.  Bull.  10, 
Sanitary  &  Eng.  Research  Project,  Univ.  Calif., 

Berke  ley  ,  1953  . 

20.  Senn,  A.,  Paleogene  of  Barbados  and  its  bearing  on 

history  and  structure  of  Antillean  Caribbean 
Region,  Bulletin  of  the  Amer.  Assoc,  of  Petroleum 
Geologists,  v  o  1 .  24,  No. 9,  September  1940. 


APPENDIX  A 


EXPLANATION  OF  SYMBOLS  USED  IN 
TABLES  V  TO  TABLES  XLII 


is  the  horizontal  distance  from,  the  shoreline  to  the  point 
in  the  sand  model  under  consideration. 

is  the  depth  to  the  interface  below  mean  sea  level  at  a 
point  X  inches  from  the  shoreline. 

is  the  height  of  the  free  water  surface  above  sea  level 
at  a  point  X  inches  from  the  shoreline. 

is  the  horizontal  distance  from  the  shoreline  of  the  point 
of  application  of  fresh  water  into  the  model, 
is  the  length  of  the  island  portion  of  the  model  measured 
at  sea  level. 

is  the  f low  rate  of  fresh  water  into  the  mode  1 . 

is  the  flow  rate  of  fresh  water  into  the  model  per  foot 

width  of  the  model. 

is  the  permeability  of  the  sand  in  the  model  under 
cons  iderat ion . 

is  the  distance  the  Y  axis  of  the  original  data  of  each 
test  has  been  transformed. 


VALUES  OF  Y  AND  Y 1  OBTAINED  VALUES  OF  Y  AND  Y1  OBTAINED 

IN  SMALL  MODEL  IN  SMALL  MODEL 
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H  =  114  inches 
L  =  122  inches 
K  =  0,7  x  10" 2  cm/ sec 
q/KH  =  0.9  x  10-3 
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LEGEND  FOR  FIGURES  14  TO  30 


mo 


Q  =  Rate  of  inflow  of  fresh  water  into  the  model 

Q2  =  Rate  of  extraction  of  fresh  water  from  model 

H  =  Distance  from  shoreline  of  point  of  inflow  of 
fresh  water. 

H  ^  =  Distance  from  shoreline  of  point  of  extraction 
of  fresh  water 

— ^  =  Stream  lines  of  flow  during  pumping 

_A _  =  Upper  diagram,  Position  of  water  table  before 

p  ump in  g . 

_  =  Upper  d  ia  gr  am ,  P  os  it  ion  of  water  table  immediately 

after  p  ump ing . 


hours 


Position  of  interface  after  the  given  elapsed 
time  since  pumping  began. 
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DISPLACEMENT  OF  SALT  WATER  WITH  FRESH  W'ATER0  RUN  NO 


INCHES 


102 


IX) 

X 

X 

u 

2: 


00 

UJ 

X 

u 

X 


LO 


UJ 

(X 

U2 

O 

U- 


VARIATION  OF  POSITION  OF  WATER  TABLE  AND  INTERFACE  WITH  TIME 


INCHES 


103 


VARIATION  OF  POSITION  OF  WATER  TABLE  AND  INTERFACE  WITH  TIME 
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DUE  TO  SIMULATED  PUMPING.  RUN  NO. 
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VARIATION  OF  POSITION  OF  WATER  TABLE  AND  INTERFACE  WITH  TIME 
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DIE  TO  SIMULATED  PUMPING.  RUN  NO. 
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DIE  TO  SIMULATED  PUMPING.  RUN  NO. 
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VARIATION  OF  POSITION  OF  WATER  TABLE  AND  INTERFACE  WITH  TIME 
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VARIATION  OF  POSITION  OF  WATER  TABLE  AND  INTERFACE  WITH  TIME 
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VARIATION  OF  POSITION  OF  WATER  TABLE  AND  INTERFACE  WITH  TIME 
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VARIATION  OF  POSITION  OF  WATER  TABLE  AND  INTERFACE  WITH  TIME 
DIE  TO  SIMULATED  PUMPING.  RUN  NO. _ U _ 
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VARIATION  OF  POSITION  OF  WATER  TABLE  AND  INTERFACE  WITH  TIME 
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VARIATION  OF  POSITION  OF  WATER  TABLE  ANI)  INTERFACE  WITH  TIME 
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VARIATION  OF  POSITION  OF  WATER  TABLE  AND  INTERFACE  WITH  TIME 
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VARIATION  OF  POSITION  OF  WATER  TABLE  ANT)  INTERFACE  WITH  TIME 


116 


VARIATION  OF  POSITION  OF  WATER  TABLE  ANT)  INTERFACE  WITH  TIME 
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FIGURE  30 

DISPLACEMENT  OF  INTERFACE  BY  RECHARGE  OF  FRESH  WATER 
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FIGURE  31 

STRUCTURAL  ELEMENTS  OF  THE 
CORAL  ROCK  AREA,  BARBADOS 

Taken  Iron.  Serin  (Ref.  16) 
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FIGURE  32 

RAINFALL  MAP  OF  BARBADOS 
Taken  from  Senn  (Ref  16) 
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FIGURE  33 

PRINCIPAL  DIRECTION  OF  FLOW  OF  STREAM  WATER. 

AND  APPROXIMATE  DEMARCATION  OF  SHEET  WATER,  BARBADOS 


Taken  from  Senn  (Ref  16) 


121 


FIGURE  34 

MAP  OF  BARBADOS  SHOWING  MAIN 
SOURCES  OF  WATER  SUPPLY 

Taken  from  Senn  (Ref.  16) 
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fir-.  A.  SCCTION  SHOWING  THC  rORMATlON  OF  SHEtT  WATER  ALONG  THE  WEST  COAST 

(no  •colo) 

WEST  •  EAST 


FIGURE  35 

FORMATION  OF  SHEET  WATER  ALONG 
WEST  COAST  OF  BARBADOS 


Taken  from  Senn  (Ref  16). 
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FIGURE  37 


VIEW  OF  BIG  MODE! 
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FIGURE  38 


VIEW  OF  SMALL  MODEL 


SECTION  THROUGH  WEST  COASTAL  GROUND  WATER  BASIN 


FIGURE  40.  TWO  DENSITY  SYSTEM  IN  EQUILIBRIUM,  SMALL  MODEL. 


FIGURE  42.  TWO  DENSITY  SYSTEM  IN  EQUILIBRIUM,  BIG  MODEL 
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APPENDIX  C 
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Determination  of  Permeabilities  of  Sands  from  Permeameter 


Tests 


A  „  Ot  t  awa  S  and 


Weight  of 

sand  in 

p  er  me  ame  ter 

=  75  9 

gr  ams . 

Int  erna 1 

d iame  ter 

of  permeameter 

<1” 

ON 

r— l 

II 

inches  . 

Length  of  sample 

=  9  .  1 

inches  . 

Head 

=  45 

inches  . 

Average  Q 

=  94 

cc/min. 

K  =  QL 
t  LA 

K  =  1.68  x  10" ^  cis/sec. 

Void  Ratio 

e  =  ^  w V  e  =  Void  ratio 

Ws  w  =  unit  weight  of  water 

_  _  qo  V  =  volume  of  soil  mass 

e  -  >j4 

Ws  =  dry  weight  of  soil  grains 


B.  Gopher 

Sand 

Weight 

of  sand  in  permeameter 

820 

g  r  ams 

Internal  diameter  of  permeameter 

=  1.94 

inch  es 

Length  of  sample 

=  8.9 

inches 

Head 

=  25.7 

inches 

Average  Q 

=  87.3 

cc/min. 

K  =  QL 
t  LA 

K  =  2.68  x  10-^  cms/sec. 
Void  Ratio 


e  =  G  wV 
Ws 


-  1 


e  =  0.48 


13d 


Determination  of  Specific  Gravity 


C  0  190 

Standard  Ottawa 

S  and 

Gopher  Sand 

+  w  +  s 

7  94 o  80 

785 . 25 

Temperature  T 

26 . 9 

2  7.0 

Wb  +  w 

701.40 

691.20 

Weight  sample  dry  &  dish 

305.63 

306.95 

Tare  dish 

156 . 6  1 

156.61 

Ws 

149 .02 

150.34 

Gs 

2.68 

2.67 

W'b  +  w  +  s  =  Weight  of 

flask  and  water 

and  s 

amp  1 e  at  T ° 

W^  -p  w  =  Weight  of 

flask  and  water 

at  T° 

(flask 

calibration  curve) 


Ws  =  Weight  of  dry  soil 

Gg  =  Specific  gravity  of  soil  particles  = 

Ws 


Ws  +  ^b+w 


Wb+w+  s 


. 


131 


Sieve  Analysis  of  Sands 

A.  C.190  Standard  Ottawa  Sand 

Weight  of  sample  =  797.7  grams 


S  ieve 

Size 

Weight  Retained  (grams) 

7a  Retained 

No  . 

30 

7  74.7 

97.4 

N  o  . 

50 

22.2 

2  .55 

No  . 

100 

0 . 2 

_  _ 

B  .  Gopher  Sand 

Weight  of  sample  =  831.0  grams 


Sieve  Size  Weight  Retained  (grams) 


N  o  . 

30 

5  .  7 

No  . 

50 

5  7  1.2 

No  . 

100 

241.2 

N  o  . 

200 

11.2 

P  an 

1 .  1 

7o  Retained 
0.68 
68.7 
29.1 
1  .35 
0  .  13 


